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CAVITATION IN LIQUID CRYOGENS 
IV - COMBINED CORRELATIONS FOR 
VENTURI, HYDROFOIL, OGIVES, AND PUMPS 

J. Hord 
I. SUMMARY 

This document is the fourth and final volume to be issued on the 
results of cavitation studies performed at the NBS. The first, second, 
and third volumes dealt with venturi, hydrofoil, and ogive experiments 
and extended the theory for correlating developed cavitation data. This 
final volume extends these experimental data and correlative expres- 
sions to the prediction of cavitating pump performance. 

It was emphasized in the previous reports that cryogenic liquids 
require less subcooling to avoid cavitation; i. e. , less Net Positive 
Suction Head (NPSH) is required for the cryogens- -liquid hydrogen and 
liquid nitrogen were used in the NBS tests. 

Developed cavity data, consisting of pressure and temperature 
measurements within fully developed hydrogen and nitrogen cavities, 
were acquired for a venturi, hydrofoil, and three scaled, quarter-caliber 
ogives. These combined data are correlated herein using the extended 
theory derived in Volume II of this report series. Maximum benefit, 
in correlation of these combined stationary-body data, is obtained by 
using the new MTWO parameter. MTWO is a liquid phase velocity ratio 
derived from two-phase flow considerations. It is significant that the 
parameters which satisfactorily correlated the individual venturi, 
hydrofoil, and ogive data also adequately correlated these combined 
data. Consequently, these correlating parameters are proven applicable 
to a variety of body geometries (two-dimensional and axisymmetric) 
that encompass internal and external cavitating flows. 



The new and old correlating expressions are extended herein to 
apply to rotating machinery. The old expressions were modified to 
account for size effects in pumps. The size- effect data obtained with 
the scaled ogives were used. Previous pump prediction work has been 
restricted to those conditions where pump geometry, flow coefficient 
and held constant while fluid, fluid temperature and 

rotative speed are permitted to vary. This study indicates that exist- 
ing correlating expressions may be used with much greater latitude- - 
only ah /AH^^ must be held constant. This is a significant extension 
of the state-of-the-art and requires development of an independent 

method for estimating K . of pumping equipment. 

c, min 

Suggestions for estimating K . are advanced herein. All 

c, min 

currently envisioned methods require detailed knowledge of the pump 

geometry. Compilation of data relating the experimental cavitating 

pressure coefficient, K . , to the calculated noncavitating pressure 

c , min 

coefficient, C^, are needed to evaluate these suggestions. 

Briefly, this study produced a simplified and precise technique 
for computing B-factors. The B-factor concept appears to be the most 
successful approach to predicting cavitation performance of pumps. 
Analysis of convective heat transfer and mass fliax limiting processes 
resulted in improved correlative — predictive expressions applicable 
to pump design and operation. Extensive developed- cavity data were 
acquired for stationary hydrodynamic bodies immersed in a flowing 
cryogenic liquid. These data were satisfactorily correlated using the 
improved semi- empirical correlative expressions. These expressions 
were extended to predict cavitating pump performance. NASA- LeRC 
liquid hydrogen pump data were used to demonstrate the adequacy of 
this extension. It was shown that the ability to predict cavitating 
performance of pumps can be significantly extended. The extension 
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of semi- empirical correlating expressions from stationary-bodies to 
pumping machinery was substantiated. Directions for future technical 
effort were discussed. 

2. INTRODUCTION 

Cavitation is usually categorized as gaseous or vaporous; in this 
study we are concerned only with vaporous cavitation. Over the past 
30 years much confusion in the cavitation literature is attributable to 
the failure of authors to clearly recognize and/or categorize their 
cavitation data. The presence of dissolved gases, minerals, and 
particulate matter adversely affects the cavitating performance of 
equipment and simultaneously complicates analysis and data correlation. 
All data treated in this report are considered to be virtually free of 
these effects- - such are the advantages of cryogenic test fluids. Vapor- 
ous cavitation is the formation of a vapor phase within the bulkstream 
of a flowing liquid, or at the interface between a flowing liquid and a 
solid surface, caused by a reduction in pressure. 

To design liquid handling equipment such as pumps and flow- 
meters, the designer must determine whether cavitation will occur 
and, in many cases, to what extent. Since the formation and collapse 
of vapor cavities alter flow patterns, cavitation may reduce the effi- 
ciency of pumping machinery [l], ^ produce vibrations, cause flow 
instabilities, and reduce the precision of flow measuring devices. 
Collapse of these vapor cavities can also cause serious erosion 
damage [2] to pump blades, flowmeter vanes, etc. We usually attempt 
to avoid cavitation but it is sometimes beneficial; therefore, it is 
desirable to learn how to predict the cavitating performance of equip- 
ment. The weight of rocket vehicles can be reduced by designing 


Numbers in brackets indicate references at the end of this report. 
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cavitating pumps to operate at higher rotative speeds and lower pro- 
pellant tank pressures; i.e., the pump is designed to swallow limited 
quantities of vapor with little degradation of performance. The potential 
decrease in pump size and weight, offered by higher rotative speeds, 
has numerous industrial and aerospace applications. 

Pumps designed to perform satisfactorily with controlled cavita- 
tion normally require an inducer upstream of the main pump rotor. 

The function of the inducer is to increase the fluid pressure enough to 
preclude cavitation in the following stage. Inducers may be an integral 
part of the pump rotor or may be separately mounted on the pump shaft. 
Inducers permit satisfactory operation of the succeeding pump machinery 
while appreciable cavitation is occurring on the suction surfaces of 
the inducer blades. Thus, net positive suction head (NPSH) require- 
ments for the pump assembly are lowered considerably with little 
sacrifice in pump performance. 

It is difficult to predict the NPSH required for the satisfactory 
performance of cavitating equipment. NPSH requirements are deter- 
mined by the complex and combined effects of fluid properties, pump 
geometry, pump operating point, and heat and mass transfer. This 
topic is treated in detail herein. The noncavitating performance of 
hydraulic equipment may be predicted from established similarity laws 
but cavitating performance is much more difficult to predict from 
fluid- to- fluid. Recent advances in this area have been made by NASA- 
LeRC personnel [3-6] and others [7-9], but additional work is required 
to improve the current technique for predicting cavitating performance 
of equipment from fluid- to- fluid. The effects of fluid properties on the 
cavitation performance of equipment are well recognized [10-19] and 
have been treated in considerable detail in previous reports [20-23]. 
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A more general predictive technique that includes the effects of equip- 
ment geometry and size in addition to fluid properties is sorely needed. 
This report outlines initial efforts in this direction and clearly shows 
that more knowledge is needed to extend our predictive capability to 
equipment design. Predicting cavitating performance for a specific 
fluid-pump combination using test data from a different fluid -pump 
combination is our ultimate goal. 

NASA has undertaken a program [1] to determine the cavitation 
characteristics of various hydrodynamic bodies and the thermodynamic 
behavior of different fluids in an effort to obtain improved design cri- 
teria to aid in the prediction of cavitating pump performance. The 
study described herein was conducted in support of this program. 

Liquid hydrogen and liquid nitrogen were chosen as test fluids 
for this study for the following reasons: (1) the ultimate goal of this 

program is to acquire sufficient knowledge to permit intelligent design 
of pumps for near-boiling liquids, and (2) predictive analyses indicated 
[1] that the physical properties of hydrogen and nitrogen make them 
particularly desirable test fluids. The objectives of this study were 
(l)to determine experimentally the flow and thermodynamic conditions 
required to induce desinent (or incipient) and developed cavitation on 
various hydrodynamic bodies, (2) to improve existing correlative 
expressions for the prediction of cavitating performance of hydraulic 
equipment, and (3)to establish, if possible, a technique for predicting 
the fluid -handling capability of different cavitating equipment using 
different fluids. The first two objectives have been effectively accom- 
plished and the third item, though ambitious, has been significantly 
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advanced. The second and third objectives are extensions of the 
state-of-the-art and more experimental pump data are needed to fully 
evaluate progress on the last objective. 

Due to the bulk of the subject matter, reader familiarity with the 
first three volumes of this report series is necessary and is assumed 
throughout this final volume. Definitions, terminology, etc., are con- 
sistent throughout this report series. Although desinent and developed 
cavity data were acquired in the stationary-body tests, only the developed 
cavity data are of interest in this final volume. 

This report summarizes the work performed on five hydrodynamic 
bodies with three different geometries (shapes). It also documents the 
analytical and correlative work performed at NBS on the experimental 
pump data [4, 6] generated by NASA-LeRC personnel. Previous cor- 
relative efforts [3-6] and those reported herein indicate that results 
obtained from stationary-body studies are directly applicable to the 
correlation and prediction of cavitating pump performance. Cavita- 
tion data for the hydrodynamic bodies (venturi, hydrofoil, and ogives), 
apparatus details, test procedures, etc., were reported in Volumes I, 

II, and III of this report series [20-22]. These three volumes also 
contained data correlation efforts for each individual body shape, 

A similarity equation, based upon the B- factor concept of Stahl 
and Stepanoff [ll], has been developed [19] for correlating cavitation 
data for a particular test item from fluid- to- fluid; this correlation is 
also useful in extending the velocity and temperature range of data for 
any given fluid. Thermal boundary layer considerations and two-phase 
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mass flux limiting concepts were used [21] to improve this correlative 
expression, A new correlating parameter, MTWO, was developed [21] 
and has proven to be a valuable correlating parameter for our venturi, 
hydrofoil [21], and ogive [22] data* The MTWO parameter is the ratio 

of V /V where V is proportional to the two-phase liquid-vapor sonic 

of. f. 

velocity across the cavity interface, see reference [21]. The correla- 
tive expressions [21] developed in the course of this study are vised to 
correlate the combined experimental data for the ogives, hydrofoil, 
and venturi. In turn, these results are used to develop cor relative /pre- 
dictive equations for pumps and are applied to the NASA-LeRC pump 
data. These generalized equations account for changes in pump size, 
geometry, flow coefficient and rotative speed, and for changes in fluid 
and fluid temperature. Prediction of the cavitating performance of 
pumping equipment was previously restricted [4-6] to changes in pump 
rotative speed, fluid, and fluid temperature. More experimental pump 
data are needed to fully evaluate the new generalized correlative/predictive 
equations. 

An improved method of predicting the cavitating performance of 
pumping equipment is presented herein. This improved method is an 
extension of previous correlative efforts [3-6], 

3. EXPERIMENTAL EQUIPMENT, INSTRUMENTATION, AND 

TEST PROCEDURES 

The experimental facility, instrumentation, error statements, 
visual and photographic aids, and test procedures used in this study 
were fully described in the first volume [20] of this report series. One 
additional error statement, concerning uncertainty in pressure measure- 
ment for nitrogen test fluid, was needed and given in section 5. 2 of 
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Volume II [21], Geometric and construction details for each body- 
t\innel configuration were given in the respective reports for the 
venturi [20], hydrofoil [21], and ogives [22], Each body-tunnel combina- 
tion was located between the supply and receiver dewars of a blowdown 
flow system as previously described [20]. 

Briefly, the objective of this experimental program was to obtain 
and analyze desinent and developed cavitation data for various hydro- 
dynamic bodies immersed in flowing cryogenic fluids. This experi- 
mental data, combined with an analytical study, was intended to advance 
existing knowledge of the cavitating performance of equipment. Informa- 
tion concerning model size, geometry, fluid, fluid velocity, and temper- 
ature is prerequisite to the improvement or development of a general 
expression for predicting the cavitation performance of hydraulic equip- 
ment. Experimental apparatus consisted of supply and receiver dewars 
equipped to control fluid state and fluid flow conditions. Instrumented 
test sections (tunnels) and streamlined bodies were situated between 
the dewars. Desinent and developed cavitation data were acquired at 
various fluid temperatures and velocities in liquid nitrogen and liquid 
hydrogen. Closed circuit television and movie film were used to define 
and analyze the data- -particularly the cavity shape data. High speed 
data recording was provided by a magnetic tape data- acquisition system. 
A short description of each hydrodynamic body is given in the following 
paragraph. 

A transparent plastic venturi [20] with a quarter- round stepped 
throat provided needed data relative to internal flow geometries. Re- 
presentative data for two-dimensional bodies were acquired for a sub- 
merged hydrofoil [21] with cylindrical (half- caliber) leading edge. 
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The three scaled, quarter- caliber , ogives [22] used in this experiment 
consisted of submerged cylindrical bodies with quarter- caliber rounded 
noses. The diameters of these cylindrical bodies varied by a factor 
of two. These ogive data supplied necessary information relative to 
size effects and external flow over axisymmetric bodies. The hydrofoil 
and ogives were sting-mounted and housed in appropriate transparent 
plastic tunnels; therefore, tests with these bodies permitted us to photo- 
graphically determine cavity shape as a function of body shape and flow 
conditions. The hydrofoil and ogive test assemblies were installed 
in the same space allocated for the plastic venturi [20] in the experi- 
mental apparatus, 

4. DATA ANALYSIS 

We seek to apply the results of this experimental and analytical 
investigation to the design and operation of liquid pumps. The NASA- 
LeRC pump impeller and inducer data [4, 6] were chosen to evaluate 
our predictive expressions. Our data and the NASA-LeRC experi- 
mental pump data are correlated and discussed separately in this 
section. 

4. 1 Selection of Stationary- Body Data 
The data selected were generated during the course of this experi- 
mental program. All of the venturi [20], hydrofoil [21] and ogive [22] 
data were used. Fluid properties for nitrogen [24] and hydrogen [25] 
are well known and documented. These property data are considered 
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more accurate than those for water [26] and common refrigerants [27] 

Other applicable published stationary-body data [3, 19, 28] were not used 
because (1) they were not available in sufficient detail or quantity, or 
else (2) they were acquired using test fluids having poorly defined fluid 
properties. The need for accurate fluid property data cannot be over- 
emphasized, as small uncertainties in property data can result in large 
uncertainties in calculated B-factor [23]--a vital constituent of our 
cor relative- predictive expr es sions . 

4.2 Correlation of the NBS Stationary- Body Data 

Previous reports [20-22] have documented correlative efforts for 
each body shape. In this final volume the combined venturi-hydrofoil- 
ogive (VHO) data are correlated. Knowledge of fluid properties, flow 
conditions, body geometry, and cavity lengths are required. The 
visual cavity lengths, as determined from movie film, were used in all 
correlative data fits. 

Pressure and temperature profiles within fully developed cavities 
were measured and are referred to herein as developed cavitation data. 

The bulkstream vapor pressure exceeds the measured cavity pressure; 
it also exceeds the saturation pressure corresponding to the measured 
cavity temperature. Therefore, the measured pressure depressions 
and the pressure depressions corresponding to the measured temperature 
depressions are called "pressure depressions." Alternatively, the 
pressure depression may be expressed in terms of its equivalent equili- 
brium "temperature depression." The state of thermodynamic equilibrium 
within the developed cavities was determined from these direct measure- 
ments of pressures and temperatures within the vaporous cavities. Cavitie 
developed on the hydrofoil [21] and ogives [22] were in stable thermo- 
dynamic equilibrium. Metastable vapor was found to exist in the central 
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and aft regions of hydrogen cavities developed inside of the venturi [20], 
Stable thermodynamic equilibrium prevailed near the leading edge of 
the cavities for all of the VHO data. Earlier analytical work [Z0,2l] 
indicated that this thermodynamic condition was sufficient for compati- 
bility of theory and experiment. 

The VHO developed cavity data are correlated [20] by using the 
experimental values of cavity pressure and temperature depressions. 
The pressure depression in the cavitated region is determined by 
subtracting the measured cavity pressure, in one case, and the satura- 
tion pressure associated with the measured cavity temperature, in the 
other case, from the vapor pressure of the liquid entering the test 
section. Profiles of measured pressure depression for all of the VHO 
data were presented and thoroughly discussed in previous reports 
[20 - 22], Pressure depression is maiximum at the leading edge of the 
cavity and zero at the trailing edge of the cavity. Maximum pressure 
depressions for the VHO data are tabulated below. 


Table 4. 1 Maximum pressure depressions for the VHO 
developed cavity data. 


Hydrodynamic 

Body 

Test 

Fluid 

Maximum Pressure 
Depression 

Venturi [20] 

hydrogen 

2 

15, 13 psi (10. 44 N/ cm ) 

Hydrofoil F21] 

hydrogen 

15, 76 psi (10, 87 N/cm^) 

Hydrofoil [21] 

nitrogen 

10. 08 psi (6, 96 N/cm ) 

Ogive [22] 

hydrogen 

13. 77 psi (9. 50 N/cm^) 

Ogive [22] 

nitrogen 

12. 89 psi (8, 89 N/cm*") 

1 — 
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The correlative expressions, developed in a previous report [21], 
were used to correlate the developed cavity data from the VHO experi- 



Equation (4-1) is an improved and extended version [21] of the simplified 
correlative expression developed by Gelder, et al. [19]. Equation (4-2) 
was derived [21] from basic heat and mass transfer considerations in 
Volume II of this report series. These expressions, along with the 
isentropic BELASH theory [23] and two least- squares data-fitting com- 
puter programs [20], were used to correlate the combined VHO data. 
These equations are used to correlate developed cavitation data in simi- 
lar test items and to predict the cavitation performance of a test item 
from fluid -to -fluid, or at another temperature of the same fluid, when 
limited test data from a single fluid are available. They are also used 
herein to correlate developed cavity data from one hydrodynamic body 
to another and to predict cavitation performance from one pump to 
another. First, we concentrate on the stationary-body data. 

Complete and detailed descriptions of the correlative technique, 
computational steps, and computer programs are given in reference [20], 
The correlative procedure, as previously described [20], can be 
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followed directly when using eq (4-1). To use eq (4-2), simply sub- 
stitute MTWO for in the computer program. Briefly, this correlative 
procedure ensures that the B values calculated from eq (4-1), or eq (4-2), 
and the BFLASH values [23] for each data point are as nearly identical 
as possible; both B values at each data point are evaluated from experi- 
mental data and this correlative procedure produces the best possible 
agreement between the isentropic flashing theory [23] and the correla- 
tive expression- - eq (4-1) or eq (4-2). This "best-fit” of the experi- 
mental data is obtained by selecting appropriate exponents for each of 
the correlative parameters in the correlative expression--eq (4-1) or 
eq (4-2). The exponent selecting process is quite complex and is treated 
in appropriate detail in reference [20]. Exponents for eq (4-1) and 
eq (4-2) were derived to evaluate the suitability of MTWO as a correlat- 
ing parameter for the combined venturi, hydrofoil, and ogive data. One 
should recall that the individual venturi, hydrofoil, and ogive data 
correlations were significantly improved [22] by use of the MTWO 
parameter; similar improvement was obtained with the combined VHO 
data. 

In eq (4-1) and eq (4-2) the cavity lengths were evaluated at the 
visually observed lengths. BFLASH was obtained for each experi- 
mental data point by using T^, the maximum measured cavity pressure 
depression, and the calculation method outlined in reference [23]. The 
fluid physical properties are evaluated at and T^ with the exception 
that MTWO is evaluated at the minimum measured cavity pressure. 

The characteristic dimension, D, was taken as the contour radius for 
each hydrodynamic body because all three body shapes possess a constant 
radius contour upon which the cavity originates. The standard deviation in 
B is computed for each set of exponents; the individual exponents may be 
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held constant or be chosen by the computer. The standard deviation in 
B factor is minimized in the computer programs v/hen one or more of 
the exponents is selected by the computer; the absolute minimum standard 
deviation is obtained when all of the exponents are selected by the computer 
as in this report. In those cases where the exponents are held constant, 
the standard deviation cannot be minimized and is merely computed. 

The set of exponents that produces minimum standard deviation in B is 
selected as the best solution for any particular batch of data; i.e., the 
standard deviation is a measure of the validity of the correlative theory 
and of the isentropic- flashing theory, as both are evaluated from experi- 
mental data* 

Because MTWO proved to be a valuable correlating parameter for 
the individual venturi, hydrofoil, and ogive data [22], the combined VHO 
data were correlated with and without the MTWO parameter to further 
evaluate its influence. Correlation of the combined VHO data with eq (4-1), 
and then with eq (4-2), provides direct evaluation of MTWO as a cor- 
relating parameter- -identical comparisons were prepared for the indi- 
vidual venturi, hydrofoil, and ogive data [22], Table 4,2 summarizes 
the individual and combined correlative results for developed cavitation 
on all geometries (bodies) tested in this study, i.e,, venturi, hydrofoil, 
and ogives. The results given in table 4.2 are discussed in the following 
section of this report, 

4.3 Discussion of the NBS Stationary- Body Data 

Cavity visual appearance, interpretation of data, and graphical 
display of typical developed cavity data were treated in detail in previous 
reports [20-22], In this final volume we concentrate on our efforts to 
correlate the combined data for the VHO. In brief, the pressure 
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Table 4. 2 Summary of correlative results for developed cavity data-- venturi, hydrofoil, ogives, and 

combined venturi-hydrofoil-ogives (VHO). 



ft Mean Percent Difference 



depressions obtained from the cavity temperature measurements were 
generally in good agreement with those derived from the measured 
pressures for the ogives and hydrofoil, i. e. , within the allowances of 
instrument error the cavity vapor was in stable thermodynamic equili- 
brium, Evidence of metastable vapor was presented in the venturi 
study [20], The VHO data revealed that cavity pressure depressions 
generally increase with increasing cavity length, body size, fluid 
temperature, and velocity; however, the pressure depressions decrease 
with increasing velocity for the two smallest ogives tested in liquid 
hydrogen. Similar results were obtained in pump inducer performance 
tests [6] using liquid hydrogen. 

These cavity parameter functional dependencies are also ahown 
by tabulating the values of the exponents derived to fit the experi- 
mental data--see table 4, 2, In reference [23] it is shown that the pres- 
sure depression increases with increasing T^ and B. Referring to 
line 4 of table 4, 2 we observe that B increases with increasing V^, 

Z, and D for the VHO data. Then, for these data, the pressure depres- 
sion must increase with increasing T , V^, t , and D, By inspecting 
the exponent data in table 4, 2, similar deductions may be drawn for the 
various body-fluid combinations tested, 

4. 3, 1 Mathematical Correlative Results 

Equations (4-1) and (4-2) were fitted with numerical exponents 
derived by forcing these equations to fit the venturi, hydrofoil, and 
ogive experimental data. These equations were derived in the course 
of this study [22] and represent extensions of the work of Gelder, et al. 
[19]. The exponents given in table 4, 2 were obtained with a least- 
squares fitting technique and a digital computer; the suitability of the 
various exponents to the experimental data is indicated by the standard 


16 


deviation in B-factor as explained previously. The mean percent dif- 
ference in B-factor, as defined in table 4. 2, may also be used to assess 
the perfection of the correlative data fit. In the hydrogen venturi data 
the value of B ranges from 2. 0 to 5. 0; in the hydrofoil experiments the 
value of B ranges from 1. 0 to 5. 0 for both hydrogen and nitrogen; in 
the ogive data the value of B ranges from 0, 7 to 2. 8 for hydrogen and 
from 0. 5 to 3. 0 for nitrogen. 

The expressions used to correlate the experimental data are 
given at the bottom of table 4, 2. The mathematical technique used to 
derive the exponents can easily pick an extraneous value for any of 
the exponents if there does not exist significant variation in the cor- 
responding physical parameter. This point has been elaborated upon 
in previous reports [20-22] and will not be discussed here; however, it 
is the lack of variation in O' that explains why El tends toward a negative 
number in line 5 of table 4, 2. In lines 1, 2, and 3 of table 4. 2 the 
numerical value of El is so small that the a term could be ignored. 

The variation in o' for the hydrogen-nitrogen correlations was about 
2:1, There was more than 400 percent change in a in the hydrogen-re- 
frigerant 114 data correlated by Moore and Ruggeri [3], and therefore the 
value for El reported in line 9 of table 4.2 is preferred when correlating 
with eq (4-1). It is apparent that the combined fluid correlations , for 
any hydrodynamic model or correlative expression, are to be preferred 
because of the greater variation in physical parameters. 

In all of our individual venturi, hydrofoil, and ogive data, use of 
the V and a terms improved the correlations; however, it is felt that 
use of these additional correlating parameters is not justified unless 
they substantially improve the correlative fit. None of the data were 
materially improved by the use of a; therefore, values for E5 are not 
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included in table 4.2. Similarly, the v term was of value only for 
some of the combined fluid correlations using eq (4-1) (see values for 
E4 in table 4. 2). Because the a term was of negligible value in the 
correlative fit, even though it varied by a factor of three, we must 
conclude that a is not an important correlating parameter for the cryo- 
gens tested; however, it may yet prove to be a valuable correlating 
parameter for other fluid combinations. Also, v may be an excellent 
correlating parameter for other fluid combinations and is of considerable 
value for some of the hydrogen- nitrogen combinations. Correlation of 
the hydrogen- refrigerant 114 data (line 9 of table 4.2) would most likely 
be improved by using one or both of these terms. 

Inspection of the first four lines of table 4.2 shows that El is 
negligible and E2, E3, and E6 are relatively consistent for the data 
correlated using eq (4-2). No such consistency exists when correlating 
the data with eq (4-1) (see lines 5 through 9 in table 4.2). That eq (4-2) 
is quite superior to eq (4-1), as a correlative expression, is readily 
shown by line-to-line comparison of the results given in table 4.2. A 
substantial reduction in standard deviation in B-factor (or mean percent 
difference in B-factor) is achieved in all cases where eq (4-2) is used-- 
compare lines 1 and 5, 2 and 6, 3 and 7, 4 and 8. The use of MTWO is 
clearly superior to the use of V as a correlating parameter. Careful 
study of the data in table 4,2 shows that the numerical values of El, 

E3, ^6 are reduced when eq (4-2) is used. The effect of the a term 

is d ecause of the many thermophysical and thermodynamic fluid 

pro embodied in the MTWO parameter. The predominant influence 

of tl lO term must also be responsible for the reductions in E3 

and E uen eq (4-2) is used. The importance of MTWO as a correlating 
parameter emphasizes that mass transfer plays an important role in 
the cavitation process. 
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Perusal of the E6 data in lines 3 and 4 reveals that the size effect 
is more pronounced for the scaled ogives than for the combined VHO 
data. Superficially, it appears that the value of E6 in line 4 should be 
favored because it correlates data for three different body shapes; 
however, each body shape has constant-radius contours. Consequently, 
it is not yet clear whether the VHO exponents (line 4) or the scaled ogive 
exponents (line 3) should be selected for future use. Both sets of 
exponents were used herein to correlate pump data and produced nearly 
identical results. The nearly negligible value of E6 given in line 9 of 
table 4.2 may be attributable to the small scale ratio (1.414:1) used in the 
refrigerant- 1 14 tests [3]. Billet [28] used water and refrigerant-113 to 
test 0, 24-inch and 0. 50 -inch diameter zero- caliber ogives. He obtained 
values of El = 0, 60, E2 = 0. 30, E3 = 0. 58, and E6 = -0. 25 for these 
tests--as derived from a formulation similar to eq (4-1). Thus, it 
appears that size effects can vary with equipment geometry. The 
experimental data represented in lines 3 and 4 are far more extensive, 
and the values of E6 are much more consistent from body-to-body than 
the earlier work [3, 28]; therefore, the size effects indicated by lines 3 
and 4 are to be preferred. 

Excluding size effects, line 9 represents the maximum variation 
in correlating parameters when eq (4-1) is used. Thus, El, E2, and 
E3 from line 9 were combined with E6 from line 8 (or line 7) to produce 
an improved expression with the form of eq (4-1) for use in pump cor- 
relation work, i. e. , El = 1. 0, E2 = 0. 8, E3 = 0. 3 and E6 = 0. 6 in 
eq (4-1). Alternatively, we could have used the exponent sets from 
lines 7 or 8 for comparative correlations of pump data. The latter 
choices were rejected because (1) the exponent sets in lines 7 and 8 
are quite inferior to those of lines 3 and 4,and(2) the El, E2, E3 
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values in line 9 have been moderately successful in previous pump 
correlations [4-6]. Equipping eq (4-1) with El, E2, and E3 from 
line 9 and with E6 from line 8 produces an improved correlative expres- 
sion of the original simplified form [4], This expression can then be 
compared with eq (4- 2)- -equipped with exponents as in lines 3 or 4 of 
table 4.2--in predicting cavitating pump performance. 

Attempts were made to improve the VHO data correlation (line 4 
of table 4.2) using new correlating parameters. Simultaneously con- 
sidering (1) the rate [20] at which a particle of liquid traverses a sub- 
saturation pressure region in cavitating equipment, and (2) the theoretical 
heat and mass transfer processes outlined in Volume II [21], several 
potential correlating parameters may be contrived. Each new parameter 
may then be evaluated by substituting them one at a time into eq (4-2) 
and correlating the VHO data. The new parameter may be inserted in 
place of the a, v, or a terms in eq (4-2). If use of the parameter 
results in a reduction of standard deviation in B- factor, the parameter 
could be retained for future use. None of the newly contrived parameters 
offered sufficient promise to warrant retention. 

V 

The parameters tried were 6, a/D, -6/ 4, K , K . , C , 

^ ' V c,min p 


V , V V 2 , '''' V 

C /K . C / 4, 4 C , C V /2g 4 , h - h , (h - h) / 4, 

p c, min , p p ’ p o c o o 

V 

^o”^^r \/Z n V V 

— 1 + (1 - C ) ' I, h - h , NPSH, (NPSH + h - h ), A , D , 

4L pjv voo 


, and (all symbols are identified in the nomenclature section 
of this report). Because the three body shapes evaluated have constant- 
radius contours the 4 /Ds« constant for all three geometries; therefore, 
substitution of 4 for D in eq (4-2) produces almost identical results. 


It is felt that 4 could be an important parameter in correlating other 
stationary-body data and in correlating pump data. In accordance with 
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the arguments presented in the dimensional analysis of Volume II [21], 
several Dimensionless Numbers were also tried as correlating parame- 
ters. None of them improved the correlative fit for the VHO data--Fr, 
We , Re , Pr, Re„» a-nd Pe were tried in eq (4-2), 

'V -{y -L' 

4. 3, 2 Influence of Body Shape, Fluid, and Size on K 

c , min 

Variations in K . from body-to-body and fluid -to -fluid were 
c,min 

previously [22] discussed in detail; however, this topic is so pertinent 

to the correlation of pump data that it will be summarized here. 

The arithmetic mean value of the developed cavitation parameter, 

K . , does not vary appreciably for the venturi data presented in 
c, min 

table 4. Z. This parameter was also relatively constant for the hydrofoil 

data [21]; this is an important result, because constant K . , eq 

c , min 

(4-1), and the isentropic flashing theory [23] are used to predict [4, 5] 

the cavitating performance of a particular piece of equipment. The 

fact that K . is different for different bodies curtails the current 
c, min 

predictive techniques [4, 5] to a particular piece of equipment; i. e. , the 
geometry (shape) of the cavitating equipment must be identical or simi- 
lar. Actually, it was anticipated [20] that min’ rnany cavitating 
bodies, would not remain constant- -as with the venturi- -for all fluids, 
cavity lengths, velocities, temperatures, etc. Then, it was neither 
surprising that min hydrofoil varied slightly [21] nor that 

Kc varied appreciably with the ogive body-fluid combinations [22]. 

For the ogive data, K . varied by a factor of 1. 5:1. 

c, min 

Also, K . varied more for the ogives and the hydrofoil than 
c , min 

for the venturi; K . was within 7 percent of K . for the venturi 
c,min c,min 

and within 15 percent for the hydrofoil. K . for the ogives varied 
^ ^ c,min ° 

as follows: 1) 0. 210-inch ogive--K . varied from 0. 33 to 0. 51 

® c , min 

with hydrogen and from 0. 36 to 0. 63 with nitrogen, 2) 0. 357-inch 
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ogive- -K . varied from 0, 50 to 0. 77 with hydrogen and from 0. 43 
— c , min ^ 

to 0, 7 1 with nitrogen, 3) 0, 420-inch ogive--K . varied from 0. 53 

c , min 

to 0. 71 with hydrogen and from 0. 45 to 0, 63 with nitrogen. From these 

data, it is apparent that niin with increasing ogive size. 

The 0. 24-inch diameter zero-caliber ogive data of Billet [28] for 

water and refrigerant-113 show similar variations in K , (0, 32 to 

c,min 

0. 52). Also, significant variations in rnin been detected in 

pump impeller [4] and pump inducer [5,6] performance tests. Then, 

for a specific piece of equipment K Yxiin vary appreciably with 

fluid and flow conditions, and K . for similar equipment can be 

c, min 

expected to vary with size. 

An unsuccessful attempt was made to establish a functional 


dependency of K . upon t , V , and T for these stationary-body 
c, min o o 

data. For the venturi [20], K . is nearly independent of t, V , 

c,min o 

and T , The hydrofoil K . is nearly independent of t and V 
o c , min / JT Q 

and increases very slightly with increasing T^ . Likewise, with the 

ogives, K . is almost independent of t and V and increases with 
c, min o 

increasing T^ , Billet obtained different results for zero-caliber 

ogives [28]--K decreased with increasing t, increased slightly 

with increasing V^, and was nearly independent of T^ . Consequently, 

the behavior of K . for different equipment is not currently predictable 
c, min 

prior to testing. 

Where K . does not vary appreciably, it is convenient to use 
c , min 

K . for predictive purposes. With large variations in K . , as 

c,min ^ i. JT o c, min 

in our ogive data, this practice will produce relatively crude predictive 
results; however, in practical applications similar flow conditions can 
sometimes be selected [4-6] so that predictions can be made at 


identical values of K . . Data presented herein shows that K . , 
c , min ^ c , min 

and consequently K varies widely with body or equipment 
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geometry, as does the pressure coefficient, C . Thus, it is quite 

P 

obvious that prediction of cavitation performance, from one piece of 
equipment to another, will require significant advances in the deter- 
mination of YYiin f^i^ctions of body geometry, fluid, body size, 

V , and T , 
o o 

As a preliminary step, we can supply, from this study and others, 

data that relates K . to the noncavitating minimum pressure coef- 
^ c, min ^ ^ 

ficient , The definitions of these two parameters are nearly identi- 

cal , except that K . is based upon minimum cavity pressure in 
^ c, min ^ ^ ^ 

V 

cavitating flow and C is based upon minimum pressure in noncavitating 
V P 

flow- -also , has a negative numerical value (see nomenclature). 
Experimental data from this study and others are plotted on figure 4. 1; 

V 

K . , K . , and C are not corrected for blockage in the upper 
c, min c,min p 

curves on this plot. In the lower curve K . has been corrected for 

c , min 

blockage and is plotted as (K . ) . If a designer can estimate values 
V c, min 00 

of C from idealized fluid flow solutions or from model scale-up tests 

P _ 

in wind tunnels, a corresponding value of K . can be picked from 

c, min 

figure 4.1. It may then be possible to apply the predictive techniques 
outlined in the later sections of this report to estimate cavitation per- 
formance from one piece of equipment to another. 

4.3.3 Developed Cavity Shapes 

One of the main objectives of the hydrofoil and ogive experiments 
was to obtain cavity volume- thickness data in an effort to improve the 
correlative theory. The hydrofoil and ogive cavity shape data were fully 
discussed in previous reports [21, 22] and are briefly summarized here. 
Cavity thickness and volume increased with increasing cavity length 
and were nearly independent of and T^. The shapes of all of the 
cavities were adequately represented over the front halves of the cavities 

p 

by a simple algebraic expression of the form 6^ = ' where 6^ and x 
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Figure 4 . 1 Variation of minimum cavitation parameters with minimum non-cavitating 
pressure coefficient, C , for various hydrodynamic bodies. 































are in millimeters and x cannot exceed The mathematical expres- 

sions [22] for the scaled-ogive cavity shapes were reasonably consistent 

0. 75 

and are approximated by 6^?^ 0.41 x . 

4.4 Selection of Pump Impeller and Inducer Data 

The NASA-LeRC [4, 6] impeller and inducer data using liquid 
hydrogen test fluid were used. Other published data [8, 9, 12, 16, 18] 
were rejected for failing to meet one or more of the selection criteria 
cited below. 

In order to correlate cavitation performance from pump-to-pump, 
we need precise information about (1) the geometry of each pump, (2) 
cavitation performance characteristics for each fluid and each pump, (3) 
instrumentation uncertainties, (4) pressure depression vs B-factor 
characteristics for each fluid, and (5) the uncertainties in fluid property 
data. Some of these selection criteria are interrelated, but each is of 
sufficient importance to warrant separate discussion. 

1) Pump geometry : Geometrical details of pumping equipment 
are usually not published along with the pump performance data. This 
may be partially attributed to protection of proprietary design informa- 
tion and partially to the complexity and burden of accurately describing 
such equipment. The pump inlet area and the impeller and inducer 
tip diameters are essential information while blade angle, blade tip 
solidity, blade configuration, etc., are highly desirable parameters 
for future correlative work. Even the highly detailed geometrical data 
published by NASA [4, 6, 30-3 2] were inadequate for some of the analysis 
performed herein- - supplementary data were obtained by personal com- 
munication [33]. Thorough documentation of pump geometry, including 
scaled drawings and photographs, will be necessary in future work. 

The NASA-LeRC publications [6,31,32] contain good examples of the 
required information. 
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2) Cavitation performance characteristics: Typical performance 

data are normally given at restricted or singular values of the flow 

coefficient 0 , To correlate cavitating pump performance from one 

pump to another requires the use of complete pump performance data« 

NPSH-AH curves at constant values of N, T , and 0 simplify the cor- 

o 

relative procedure. From these data it is possible to determine 
and then for a prescribed value of i|j'. Low values of (approxi- 

mately 0. 7) may be picked to provide the best available NPSH data. 

Higher values of \|;^ require more precise measurement of NPSH and 
more precise determination of AH^^ because the slopes of the NPSH - AH 
curves are flatter at higher values of ^--typical NPSH- AH curves are 
shown on figure 4, 2, Most published data use very high values of 
(0. 97 or 0. 98) which effectively amplify the uncertainties in measured 

NPSH and AH . Complete NPSH - AH data are available for the 
nc 

NASA- LeRC pump impellers £4,31,32] and inducers [5, 6, 34-36], 

3) Instrumentation uncertainties : The precision of pump perform- 
ance data is frequently difficult to assess because instrument uncer- 
tainties are not clearly stated. Knowledge of the uncertainties in the 
measurement of NPSH, AH, T , N, and flow is required to perform 
accurate correlations or predictions. Performance data with cryogens 
or conventional liquids, such as water, require very precise measure- 
ment of temperatures and pressures. The steep slope of the saturation 
pressure-temperature curve for water signifies that small errors in 
temperature measurement can introduce large errors in NPSH, Pre- 
dictive work with water is further complicated because cold water 
exhibits low values of Ah ,and vaporous cavitation effects are some- 

V 

times masked by gaseous cavitation. Thus, instrumentation require- 
ments [9] for water correlations are very exacting. It is common 
knowledge that precise thermometry instruments are needed when cryo- 
genic test fluids are used. 
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4) Ah vs B-factor: The B-Ah characteristics of each fluid 

V V 

are markedly different [23], varies almost linearly with B-factor 

for water, but Ah is a nonlinear function of B for the common cryo- 

V 

gens [23], Cold water and other conventional liquids, with boiling 
points near room temperature, characteristically exhibit relatively 
low values of Ah . In contrast, cryogenic fluids show relatively large 

V 

values of Ah . Fluid properties must be accurately known to produce 
accurate B - Ah data, 

V 

5) Fluid property uncertainties: Fluid property data are simply 

inadequate for most conventional fluids --particularly hydrocarbons 
2 

[9, 16, 18] ; water properties are in many respects still marginally 

defined even though much effort has been expended in this area [26, 

37, 38], A discontinuity in the saturation temperature-entropy data 
for refrigerant- 1 1 3 [27] does not encourage its use as a test fluid. As 
a result of stringent aerospace and cryogenic industry demands, thermo- 
physical properties of most pure cryogenic fluids [24, Z5, 39-43] are 
better known. Because of particle-physics and aerospace emphasis, 
the transport and thermodynamic properties of parahydrogen [43,44] 
are probably the most extensive and precise. It was previously shown 
[23] that small uncertainties in fluid properties can result in very large 
uncertainties in the calculated B-factor. Large uncertainties in B-factor 
result in correspondingly large uncertainties in Ah and result in crude 
correlative and predictive analyses. 


See communications by A. B. Bailey (p. 57), F. C. Gilman (p. 58) 
and M. Sutton (p. 63) at end of Chivers [9] paper. 
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If we were not attempting to correlate data from pump- to-pump, 
we could relax our selection criteria. Correlating data for a given 
pump and fluid at various inlet conditions or for a given pump with 
different fluids is much less difficult. In the latter two cases, the first 
selection criterion is not applicable and the second and third criteria 
are relaxed somewhat because a single test apparatus is used, i.e., 
the test procedures, data uncertainties, data analysis techniques, experi- 
menter judgments, etc., should be reasonably consistent. The fourth 
and fifth selection criteria are of vital importance in all correlative work. 
Selected data from Saleman [l6] and Spraker [18] were previously corre- 
lated [5] using the simplified expressions reproduced herein. These 
data were only suitable for fluid- to- fluid correlations and NPSH re- 
quirements were generally predicted within the estimated uncertainties 
of the experimental NPSH data. The predictions are undoubtedly well 
within the combined uncertainties in pump data and fluid properties data. 

After careful study of the NASA-LeRC pump data [4-6, 34-36] 
and discussion of these data with one of the NASA authors [33] , addi- 
tional selection criteria were imposed, 

a) Some of the low 0 inducer data were rejected because of the 
'dips' and 'wiggles' in the ^ - NPSH curves ( see figure 4, Z). Such 

data complicate the determination of AH and the NPSH for pre- 

nc 

scribed values of Because of the difficulty in establishing ^ 

nc 

for some of the higher 0 data, the higher values of were not used; 
i. e. , the uncertainty in graphically evaluating NPSH for a specific 
is minimized at lower values of \(r^--again, refer to figure 4. 2. All 
inducer data were restricted to values of = 0. 7. Values of i|;' = 0. 7, 

0. 8, and 0. 9 were used for the impeller data. 
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b) All NPSH data were selected so that there was minimal possi- 
bility of the existence of vapor in the pump inlet liquid, see figure 4. 2. 

c) In the data reported by NASA-LeRC , it is implied that K 

c, mm 

does not vary with T or N, but is a function of 0 only at constant 

o 

values of In pump-to-pump correlations it is convenient for K 

^ ^ ^ ^ ^ c, min 

to be a function of but two parameters, 0 and \|;^ Thus, it was imposed 

that K . p- constant (± 20 percent) for fixed values of 0 and at 
c, min 

all values of T and N for each impeller and each inducer. All of the 
o 

NASA pump data met this criterion except for the 78° inducer data 
[361, The 78° inducer data had variations in calculated K 

c , min 

that varied by 300 to 400 percent. It is felt that this selection criterion 
reflects on the internal consistency of the cavitation data for each piece 
of pumping equipment; therefore, the performance data for the 78° 
inducer were rejected. The apparent inconsistency in these 78° inducer 
data may possibly be attributable to the low cascade solidity (1. 856) 
of this inducer, A solidity ^ 2. 0 to 2. 5 improves suction performance 
and tends to counteract cavitation-induced oscillations -- see Jakobsen [45], 
The pump impeller and inducer data selected from the NASA 
publications [4-6, 34-36] are listed in table 4, 3. These data are cor- 
related herein. Geometrical details of the NASA pump equipment are 
summarized in table 4. 4 and photographs of the NASA pump impellers 
and inducers are shown on figure 4, 3. 

4. 5 Correlation of Pump Impeller and Inducer Data 
After fully evaluating the NBS stationary-body data, it was 
decided that the following two equations should be used to correlate 
the NASA-LeRC hydrogen pump data: 
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Table 4. 3 Selected data from NASA-LeRC hydrogen pump experiments 
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Denotes reference data used to compute the results recorded in Tables 4.6 and 4.7. 







Table 4. 3 Selected data from NASA-LeRC hydrogen pump experiments, (Concluded) 
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Denotes reference data used to compute the results recorded in Tables 4.6 and 4.7. 
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main blade hub -to -tip ratio 
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To show true size effect [22] , the I and D terms may be combined to 
form the dimensionless ratio, ■t/D, as follows; 



( 4 - 5 ) 


B 


B 


ref 


MTWO \ 

MTWO , 1 
ref / 


V- 

‘"'/“'refj 



( 4 - 6 ) 


The T/D ratio is used in a wide variety of geometric scaling problems, 

and the true ’size effect' is indicated by the value of the exponent on the 

D term; i. e. , for the same l/D ratio the size effect is mathematically 

0 9 0 71 

represented by * or * in eqs (4-5) and (4-6). 

The first three exponentiated terms of eq (4-5) are those used 
by Moore and Ruggeri [3-6]--the D term was derived from the NBS 
stationary-body data. Thus, eq (4-5) is used herein to provide conti- 
nuity with the cor relative /predictive work of Moore and Ruggeri and to 
provide a comparison of eq (4-5) and eq (4-6) when used to correlate 
pump data. Equation (4-6) was far superior to the mathematical form 
given in eq (4-5) for the correlation of stationary-body data; therefore, 
eq (4-6) was expected to excel when correlating pump data with these 
two equations. The results given herein confirmed this reasoning. 

To apply eq (4-5) or eq (4-6) to the correlation of pump data, 

we must relate B to Ah and then relate Ah to the NPSH requirements 

v V 

of each pump. Thus, an additional mathematical relationship is 
required and obtained from definition of cavitation parameters. First 
the cavitation parameter is defined [4] in terms of the minimum cavity 
pressure by the relation 
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(4-7) 


K 


c, min 


h - h 
o c , min 

7 - i4 + 

V /2g 

o c 


Ah 


V ^/2g 

o 


where K is the conventional cavitation parameter defined bv 

V ^ 


K 


h -- h 

o V 


^ V ^/2g 
o c 


(4^8) 


It is apparent from eq (4-7) that K . = K when Ah =0 and by defini- 

^ c,min V V 

tion we have replaced h in eq (4-8) with h . in eq (4-7). Next, 

V c,min 

we take the definition of 


V 

NPSH = h -I- - h . 

o 2g V 


(4-9) 


Combining eq (4-8) and eq (4-9) , we get 

'2g 


K 


NPSH 


V 


J 


1 


(4-10) 


Then eq (4-7) and eq (4-10) may be combined to lorm 

^ 2 

A -f K . ^ = I NPSH f Ah \ . 

V C,miny2g^ vy 


(4-11) 


Equation (4-11) indicates the relationship that must exist between NPSH 

and the cavitating pressure coefficient, K . ; it also provides the 

c, min 

framework for the second required mathematical expression. We may 
rewrite eq (4-11) in terms of pump parameters as 

.. 2 

- f NPSH + Ah^ j , (4-12) 


/ \ 

y ' ^c,minj 2 
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where = 0ND^ , Expressing eq (4-12) in,.terms of reference data, 
we obtain the desired expression 


(I +K . \ / 0ND 

_V c, min / / t 

(■ +K^,min)ref \(«"°t)re£ 


(NPSH + Ah^^ 

/^PSH 4 Ah \ ' • 

V v/ ref 


Equation (4-13) may now be paired with eq (4-5) or eq (4-6) to 
form the equation pairs needed to correlate the pump data. Equation 
(4-5) and eq (4-13) constitute our first equation pair 

r S 



(4-14)J 


where D 

t 

V = 0ND 
o 


replaces the characteristic dimension D in eq (4-5) and 
t “ 


Equations (4-14) differ from those formulated by Ruggeri and 
Moore [3-6] only in the value of the exponent on the term in the 
second equation. To assure dynamic and geometric similarity of 
cavitating flow, eqs (4-14) were restricted [3-6] to a single pump. 
Physical reasoning [4] and venturi tests [46] indicated that nain 
could be expected to vary with 0 for any specific pump. Subsequent 
pump experiments [4-6] substantiated this idea. Variation in 0 causes 
a change in flow incidence angle. Therefore, variation in 0 is analogous 
to an actual change in blade geometry for a specific pump and results 
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in a different K .for each value of 0 . Because K . was not 
c,min c, min 

predictable prior to testing, Moore and Ruggeri simplified eqs (4-14) 

by specifying that constant values of K . , 0 , and \lr^ be used to 

c , min 

predict cavitating performance of a given pump. Constant 'i/' relates 

to cavity volume and constant values of as will be discussed later 

in this report. Holding K . , 0, t(t^and pump geometry constant, 

c , min 

we obtain the simplified equation pair used with good success by Moore 
and Ruggeri [3-6]: 



B 


B 


ref 


^NPSH + Ah^) 
^NPSH + 




The third equation pair is obtained by using eq (4-6) and 
eq (4-13): 


(4-15) 
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t, ref 


(4-16) 


Referring to eqs (4-16), we note that the first equation is 
derived from basic definitions and is not restricted to a singular 
pump geometry; therefore, the first equation of (4-16) may be extended 
to predict cavitating performance from one pump to another if all of 
the parameters in this equation are known. The next section of this 
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report elaborates on this topic. For the second equation in (4-16) to be 
applicable to varying pump geometries, we must satisfy ourselves that 
the basic cavity shapes and heat and mass transfer phenomena do not 
vary markedly from one pump to another. The cavity shape data obtained 
in our stationary-body tests [22] and those from Rouse and McNown [29] 
indicate that cavities formed on a wide variety of hydrodynamic bodies are 
generally parabolic. Also, the exponent data for the second equation in 
(4-16) are remarkably consistent for the various stationary-bodies tested 
in this study (see table 4, 2 ) . This consistency of exponent values indicates 
that the heat and mass transfer processes associated with cavitation [21] 
do not vary appreciably for these hydrodynamic bodies. Thus, it appears 
that the second equation of (4-16) is equally applicable to varying pump 
geometries provided that this equation is initially evaluated from experi- 
mental data where there exists significant variation in body geometries. 

It is felt that this latter criterion is met by using the results of the VHO 
data correlated herein. 

Arguments similar to those presented above may be advanced in 
favor of using eqs (4-14) to correlate cavitating performance from pump 
to pump; however, the exponent data for the second equation of (4-14) 
are relatively inconsistent for different body geometries (see table 4, 2), 
and these arguments are less convincing. For a specific pump geometry, 
equation pairs (4-14) and (4-16) are equally applicable. Equations (4-14) 
are applied herein (to specific pumps and to varying pump geometries) to 
provide a comparison of proven predictive expressions with the improved 
expressions given in eqs (4-16), 

Equations (4-14) and (4-16) were used to correlate the NASA-LeRC 
pump data. The correlative results are summarized in tables 4. 5 through 
4. 7 and are discussed in detail in the following section of this report. 

Exact computational procedures with numerical examples are illustrated 
for eqs (4-16) in appendix A. The detailed illustrations [4] for use of 
eqs (4-15) mav also be helpful to the reader. 
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These tedious and sometimes arduous computations were computer-coded and 
supplied by William R. Parrish. 
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These tedious and sometimes arduous computations were computer-coded and 
supplied by William R, Parrish. 

















For pump designers, it is instructive to demonstrate the relation- 
ship between K . and the cavitating suction specific speed S . 

c , min c 

First, we note that the noncavitating suction specific speed 


S = (NQ^^^) / (NPSH)^'^^ 
nc 


(4-17) 


By replacing NPSH in eq (4-17) with (NPSH + Ah ) as we did in making 

the transition from K in eq (4-10) to K .in eq (4-11), we obtain 

V c,min 


S = (NQ^"^^) / (NPSH + Ah ) ^^^ 

C V 


(4-18) 


Combining eq (4-11) and (4-18) and substituting Q = A^V^, = 0ND^, 

and A = (ttD ^) /4 yields 


[(n/4)^^" (Dp/Dt) 


0 /I + K . ) 
V c,min/ 


3/4 


(4-19) 


The significance of holding 0 , K . and pump geometry constant 

c , min 

as in the development of eqs (4-15) is now apparent--S is also held 

c 

constant. Holding S = constant implies similar suction performance 
which is precisely the optimum condition for application of eqs (4-14) 
and (4-16). 

If preferred, we may now use eq (4-19) to rewrite eq (4-13) in 
terms of S^. Eliminating rnin combining these two equations 
produces 


V)ref 
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Equation (4-20) may now be paired with eq (4-5) or eq (4-6) to form 

the desired correlative/predictive expressions. These equation pairs 

will produce mathematical results that are identical to those obtained 

from eqs (4-14) and (4-l6), respectively. Again, a method of evaluating 

S (or K ) is needed in order to apply these predictive expressions 

c c, min 

to varying pump geometries. 

4,6 Discussion of Pump Impeller and Inducer Data 

Correlation of the NASA-LeRC pump data with eqs (4-14) and 

(4-16) requires K . vs 0 data for each pump. K . was calculated 
c; min ^ c, min 

for each value of 0 and for each pump using eq (4-14) or eq (4-l6) and 

the NASA-LeRC experimental data. The procedure for calculating 

K . and corresponding values of S from experimental data is 
c, min c 

described in appendix A. Calculated values of K . and S are 

c, min c 

given in table 4.3. In practice, K . (or S ) will be unknown prior 

c , min c 

to pump testing; therefore, use of eqs (4-14) or (4-l6) for predictive 

purposes will require an independent method of estimating K . . 

c, min 

Potential techniques for determining K . (or S ) are suggested in 

c,min c 

this report. 

Cavity lengths may be calculated if cavity vapor volumes and 
characteristic dimensions of each pump are known; however, we have 
insufficient knowledge concerning the extent of cavity vapor volumes 
within pump blade passages and must take - constant in eqs (4-14) 

and (4-16). This constraint is defensible only for similar cavity shapes; 
therefore, it is considered valid for homologous pumps but perhaps 
inaccurate for dissimilar pumps. As a result of increased vapor block- 
age, pump head will decrease as cavities grow within the blade passages; 
thus i' must be functionally related to for all pumps, decreases as 

V 

increases and these trends lead us to suggest a simple functional 
relationship 


cc 1 - 

V 


(4-21) 


Because is mathematically related to I/D we can also assert that 
V’ t 
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^ is a function of t/D . Obviously these functions, ^ - f ^ , 

t V 

^ ~ f (<t/D ) and ilf ^ f (f./D ), may differ from one pump to another; 

V t t 

however, holding \[r ^ = constant seems sufficient to assure that 

- constant for similar cavity shapes in similar pumps. We will 
also have to be content with this constraint for pump-to-pump cor- 
relations, even though we are fully aware that it may be inappropriate. 
Fortunately, is ^ weakly exponentiated term in eqs (4-14) and 

(4-l6) and thereby forgives many weaknesses of the imposed constraint. 

The results given in table 4. 5 are for the simplified and optimum 

condition where ifr' 0, and K . are invariant. These criteria also 
^ c,min 

require that the pump geometry be invariant because each pump was 
tested over a different range of 0’s. Thus, the reference data and 
predictions relate to identical pump geometries as witnessed by the 
numerical entries along the horizontal lines of table 4. 5, These data 
indicate the suitability of eqs (4-14) and (4-16) for predicting NPSH under 
the prescribed conditions. All possible pairs of experimental data points 
were used to establish reference data for each pump--two experimental 
test points are required to establish reference data (see appendix A)-- 
and then each set of reference data was used to calculate the mean 
percent error (MPE) in predicted NPSH. The MPE values reported in 
table 4. 5 are the arithmetic average of all possible permutations of the 
available reference data. Using a single set of reference data produces 
results similar to those given in table 4. 5 but with less statistical signifi- 
cance. Due to the scarcity of constant 0 data it was felt that the results 
shown in table 4. 5 were most representative of the available data. 

The results given in table 4.5 indicate that eqs (4-14) and (4-l6) 
are quite adequate for correlating the impeller data but that eq (4-1 6) excels 
in correlating the inducer data. The overall correlative-predictive re- 
sults are given in the bottom right-hand corner of this table. The overall 
success of the predictive expressions is assessed by calculating the 
overall mean percent error (OMPE) in predicted NPSH, 
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OMPE 


X!I(mpe) ntpp: 

^(NTPP) 


(4-22) 


Under the prescribed conditions, we could reasonably expect to predict 
NPSH within about ± 8 percent using eqs (4-16), 

Table 4. 6 lists the results for pump-to-pump predictions with 
i|r' constant; here, appropriate blocks in the table identify the reference 
and predicted p\imp geometries. As an example, the use of eqs (4-1 6) 
and reference data from impeller A at ^ = 0.7 produces MPE values of 
13.73, 5.34, 9.55, and 10.70 in predicting the cavitating performance of 
impeller A, impeller B, inducer A, and inducer B, respectively. The 
correlative equations producing the lowest MPE are listed first in each data 
block within table 4.6. The values of MPE shown in this table were cal- 
culated using the reference data identified in table 4.3. These reference 
data minimized the OMPE for eqs (4-16) and produced nearly optimum 
results for eqs (4-14). Thus, the predictive results shown in table 4.6 
may be considered the best possible using eqs (4-16), and under the most 
favorable conditions we might predict NPSH within about ±12 percent 
(see OMPE for eqs (4-16) in table 4.6). That eqs (4-l6) are clearly 
superior was again demonstrated by also computing the OMPE for all 
possible permutations of the available reference data. The resultant 
values of OMPE were 19.92 and 24.44 for eqs (4-16) and (4-14), respec- 
tively, Under the conditions prescribed in table 4.6, but using all avail- 
able reference data, we could expect to predict NPSH within about ± 20 
percent using eqs (4-l6). 

It is interesting that the values of OMPE obtained using eqs (4-14) 
do not vary much as correlative restraints are relaxed in going from 
table 4. 5 to table 4, 6 (using cyclic permutation of reference data the 
OMPE varies from 23. 09 to 24.44). Conversely, the OMPE for eqs 
(4-16) is much lower when pump geometry is held constant as in table 4. 5 
(using cyclic permutation of reference data the OMPE varies from 7. 56 
to 19. 92). Direct comparison of the OMPE values in table 4. 5 with the 
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optimum values in table 4.6 is also interesting- -the trend exhibited by 
eqs (4-l6) seems the most logical. 

In preparing the results shown in tables 4, 5 and 4. 6, it became 
apparent that eqs (4-16) dominantly produce the lowest and most con- 
sistent values of OMPE for each set of reference data. As the ref- 
erence data were varied, the individual values of MPE were also most 
consistent when eqs (4-16) were used. This improved consistency trend 
is attributed to use of the MTWO parameter in eqs (4-16), Similar 
results were obtained for the stationary-body data (see table 4.Z), 

These observations, exemplified by the results given in tables 4. 5 and 
4.6, clearly demonstrate that eqs (4-1 6) should be used for all future 
correlative and predictive work. 

As expected we obtain the best correlative and predictive results 

when if \ 0, K . and pump geometry are invariant; however, the NPSH 
c, min 

predictions are surprisingly good when only iji is held constant. Care- 
ful study of the data in tables 4. 5 and 4.6 reveals that there is room for 
improvement in the formulation of the predictive expressions. In an 
effort to improve our predictive capabilities, several ideas were explored 
as explained in the following paragraphs. 

Pump designs [14,45] employ blade relative velocities to good 
advantage; therefore, w^ was substituted for in formulating eqs (4-l6)-- 
no correlative improvement was gained. Likewise, no improvement 
resulted from substitution of blade tip chord length (C^) for in eqs 
(4-16). Substitution of the modified cascade tip solidity (C^^) for T/D^ 
in eqs (4-l6) also failed to improve the NPSH predictions. Cavitation ori- 
ginates at the tips of the main and splitter blades of the impellers. These 
blades have different radial and axial coordinates, and selection of appro- 
priate characteristic dimensions is complicated. The main blade tip dia- 
meters were selected and used herein, but it is recognized that the primary 
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and secondary splitter blades also contribute to the overall cavitation 
performance of equipment. 

Accordingly, other characteristic impeller and inducer dimensions 

were formulated and substituted, in turn, for in eqs (4-14) and (4-16). 

D , D D and D ^ were calculated (see table 4. 8 ) and sequen- 
D ml mZ m3 

tially substituted for D^: All of these D's were superior to in the 
impeller-to-impeller correlations. p^^oduced the 

lowest values of MPE as recorded in table 4. 7-- identical reference 
data are used in tables 4. 6 and 4. 7 so that the results shown in these 
two tables are directly comparable, was generally superior in all 

other correlations where the inducer data were involved, i. e. , in inducer- 
to-impeller , induce r- to -induce r , or impelle r- to- indue e r correlations. 
Because the latter groups comprise the bulk of the pump data, pro- 
vided the best overall predictions; i. e., produced the lowest values 
of OMPE. For impeller-to-impeller correlations the future use of 
D ^ is suggested. Better correlative results might possibly be obtained 
by selection of some other characteristic impeller dimension. 

Inducer suction performance is characterized [45] by the ratio 
6^/P^ where 9^ = (tan % - p/Z) and P is the blade angle. As shown 

in section 5. Z of this report ,9^ /P can also be considered an important 

. -1 

parameter in the analysis of blade cavitation- - 9^ = (tan 0 - P). For 
small values of 0 and P these two parameters, 0^/P and ^ y are 
analogous to the normalized 0 parameter (0/tan (p/Z))used by Moore 
and Meng [6]. Attempts to correlate the pump impeller and inducer 
data at constant values of 0 /p and/or 0 /P were disappointing- -bette r 
results were obtained by unrestricted variation of 0 as recorded in 
table 4. 6. 
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Efforts to improve the pump correlations by using different sets 

of exponents from table 4. 2 were also mostly unsuccessful. Equations 

(4-14) were equipped with the VHO exponent data from line 8 of table 4, 2 

and eqs (4-16) were evaluated using the ogive exponent data from line 3 

of table 4. 2, The latter approach is only slightly inferior to the use of 

eqs (4-16) as formulated, i, e. , when equipped with the VHO exponents 

from line 4 of table 4, 2. This result was anticipated because the 

exponent data for the ogives and VHO were reasonably consistent. 

Under the conditions specified in table 4. 5, NPSH predictions were 

improved by equipping eqs (4-14) with the exponent data from line 8 

of table 4. 2. The use of (V^) * etc. , with eqs (4-14) is equivalent 

to stating that ^orB?^ constant when , 0, K . and pump 

^ ref ^ c,min ^ ^ 

geometry are held constant. Equations (4-14) with B = constant produced 
an OMPE of about ± 10 percent under table 4, 5 conditions; however, the 
OMPE under table 4, 6 conditions was drastically increased by holding 
B constant. 

These findings are significant because it appears that the simplified 
predictive techniques- -cor relative restraints specified in table 4. 5 and 
use of eqs (4- 14)--proposed by Moore and Ruggeri [3-5] can be improved 
by merely using the definition of niin holding B constant, 



f NPSH + Ah ) 

^ v^ I 

(NPSH + Ah \ . 

' V' ref / 


B _ = constant 
ref 


(4-23) 


Much better results are obtained in all cases by the use of eqs (4-16). 
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5. POTENTIAL METHODS OF PREDICTING K 

c, min 

AH of the pump NPSH predictions presented in table 4, 6 are 

based on the premise that K . for S ) can be estimated without 

c , min ' c 

the benefit of pump test data for each specific pump. This section 
suggests potential techniques for estimating K , min. 

V 

5. 1 Estimation of K . from C Data 

c , min p 

(Pump experiments are not required) 

Recognizing that K vary with and may vary with 0, 

we will initially settle for rough estimates of K . ; therefore, values 
c, min ^ 

of K .at mean values of \lr^ are sought. If the ^ . vs C 

c, min ^ c, min p 

data for our stationary-bodies are valid for pumping machinery, we 

can estimate K . from figure 4. 1, This value of K . may 
c, min ® c, min 

correspond to some finite but currently unknown mean value of 

V 

To use the data on figure 4. 1, we must first determine C for each 

P 

specific impeller or inducer blade. The latter is a complex task and 
presents a challenge to the fluid dynamicist; however, a body of informa- 
tion [47-51] exists for these ideal fluid flow calculations. For very 

v 

small flow incidence angles, values of may be estimated from any 
convenient compilation of air foil data [5Z, 53] , Larger flow incidence 

V 

angles and/or the requirement for more precise values of C require 
the use of computerized numerical solutions to account for the variation 

V V 

in C with changing 0 . Once C is obtained, figure 4, 1 can be used to 
P p 

get K . . Considerable work is required to test the validity of this 
c , min ^ ' 


approach, but it appears that the effort would be justified. 

Another simple approach for estimating xnin deduced 

from a rule of thumb suggested by Jakobsen [45], He suggests that the 
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conventional cavitation number, K , can be estimated from 

V 


1 + K = 

V 


V 

C 


(5-1) 


To maintain consistency with our predictive formulations, we automa- 
tically replace K with K . and there results 
^ V c , min 


1 + K 


c, min 


V 

C 


(5-2) 


This simple rule of thumb works exceptionally well for the NBS 

V 

stationary-body data. Using eq (5-2) and the values of C tabulated 

P — 

on figure 4, 1, the reader can readily verify that values of K 

® ^ ^ ^ c, min 

are approximated rather well. It is not at all surprising that eq (5-2) 

V 

fails to reproduce Billet s data as recorded on figure 4.1. The for 

Billet' s zero-caliber ogive data is an experimental value taken from 

V 

Rouse and McNown [291. The theoretical value of C for a zero- 

L p 

caliber ogive is infinite; therefore, any deviation from a true zero- 

V 

caliber will lower the value of C . Such deviations may result from 

P 

minute cavity formations or manufacturing imprecisions. Of course, 
real fluids are also incapable of infinite accelerations around sharp 

V 

corners. Thus, the true value of C for zero-caliber ogives may well 

P 

be higher than 0. 64, It may be wiser to use eq (5-2) and the experi- 
v 

mental value of K . from Billet to estimate C for his blunt-body 
c,min p V 

cylinders. Performing this algebra, we obtain ^ - 1,40 for 

r 


3 ^ 

K and C may also be written f45l in terms of w rather than V 
V p t o 
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4 V 

zero-caliber ogives-- a tenable value when extrapolating the C data 

P 

of Rouse and McNown [29]. 

Again, the validity of eq (5-2) should be tested by application to 

V 

pumping machinery. Once C is obtained for each value of 0, eq (5-2) 

P 

may supply sufficiently accurate values of K . --again, at a yet to 

c, min 

to be determined mean value of 

V 

5.2 Estimation of K . from C Data (Compilations of ex- 

c, min p ^ 

perimental min similar pumps required) 


If the foregoing approaches fail to produce reliable values of 


K .for pumping equipment, it may be necessary to compile K 


c, min 


V 


c, min 


vs 0 and C vs 0 data for various classes (families) of pumps. Such 
P 

compilations of experimental and theoretical data would permit the 

V 

development of functional relationships between K . and C 

c,min p 

for pumping equipment. The following paragraphs outline the type of 
data needed. 

Moore and Meng [6] plotted K . vs 0 for inducers A, B, and 
^ c,min 

C--these data are shown on figure 5. la. These curves do not reveal 

any obvious trends with inducer geometry. To account for blade 

fairing angle, 0 was normalized and the same data replotted [6]- -see 

figure 5. lb. Still, no geometric trends are apparent. If we replace 

K . with (K . ) , we obtain the results shown on figure 5. 2. 

c,min ' c,min'w ^ 

The original data now resemble a family of curves with definite geo- 
metric dependencies. Improvements in the K . vs 0 data for 

c , min 

impeller A can also be demonstrated by using (K . ^ • This 

' c,min'w 

observation merely suggests that (K ^ ) may be a good parameter 

to use in compiling future experimental data. If (K . ) is used, 

y y ^ * min w 

we should also replace C with (C \ 

p ' p' w 


See p. 25 of reference [29]. 
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<#> 


o) Flow Coefficient 



Figure 5. 1 Comparison f6] of values for three helical inducers (ijf' = 0. 7). 
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( Kc,min)w ^ ^c, min) 


3t 




Figure 5. 2 Comparison of (K . ) values for three helical inducers (ijr' = 0. 7). 
® c, rmn w 
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Figure 5. 3 illustrates a typical inducer blade with the inlet 
velocity triangle constructed at the blade tip. We shall refer to this 


V 


figure frequently in that which follows, and ) is defined to apply 
only to the suction surface of the leading edge of the impeller or inducer 

V 

blade at the tip. characterizes each specific piece of equipment 

and is independent of the inlet pressure, temperature, velocity, and 

V 

the fluid. The (C j of a hydrodynamic body is constant when the fluid 
V p/w ^ 

flow is sufficiently turbulent (Re > 10 ) and the flow incidence angle is 
constant. In pumping equipment Re is normally high, but the angle of 

V 

flow incidence varies with the flow coefficient 0 . Thus, the \ 

' p' w 

for a specific piece of equipment may be considered to vary only with 0. 

V 

(Vw must be invariant if we expect identical cavitation charac- 
teristics, so for a prescribed geometry and it is sufficient to hold 
0 constant. In correlating from one piece of equipment to another, the 
(C ) profiles for the two bodies should at least be similar, even though 

^ FW V Y . 

the (C ) values may be different. The ) values may differ 

' p' w ' p' w 

because of different tip contour, fairing angle, or cascade solidity at the 
tip. To obtain similar (C^) ^ profiles , we must hold 0^ or 0^ constant-- 
either angle may be used depending upon our choice of the blade suction 
surface or pressure surface as a datum plane. 

Because the leading edges of most impellers and inducers are 
rounded aind then faired, we may consider all of them as wedge-shaped 
bodies with cylindrical noses. While each wedge angle (fairing angle) 


V 


would produce a different ) , the (C ) profiles would be quite 

'p'w 'p'w 

similar for the small blade angles under consideration and for constant 
0 . It appears that 0 is a most pertinent correlating parameter that 

accounts for both flow coefficient (0) and blade angle (3). Instead of 
holding 0 constant for predictive purposes, we could hold 0^ constant. 
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Holding 6 constant prescribes specific, but not identical, values of 

Ld 

0 for each pump. Attempts to correlate the NASA-LeRC pump data 
holding 0^ constant were disappointing. Better results were obtained 
with no restrictions on 0 as reported in table 4. 6. 

Referring to figure 5, 3, let us consider rotation of a constant 
velocity vector, w , around the blade tip. Note that w sin 0 is 
normal to the pressure surface datum plane and w cos 0 is parallel 
to this datum. Clockwise rotation of w : As 0 increases the liquid 
impacts on the suction surface, reducing the peripheral velocity 
around the leading edge of the blade. This produces higher pressures 

V V 

on the suction surface: therefore, h - h and (C ) are lowered by 

o ' p'w 

clockwise increase in 0^ , Counterclockwise rotation of w^ : As we 

increase 0 in a counterclockwise direction from the pressure surface 
datum, the liquid must accelerate to flow around the leading edge of 

V 

the blade. Pressure on the suction surface is lowered causing h - h 

V ^ 


and 




w 


of the suction surface to increase. 


Holding w^ constant oLnd rotating it about the pressure surface 
datum plcine requires a continuous change in 0; clockwise rotation of 
w^ results in larger values of 0 and counterclockwise rotation produces 
smaller values of 0 . The location of the minimum pressure point on 


the suction surface, and the value of 


<Vw 


for the suction surface, will 


change with each value of 0 for a particular pump. As 0 increases- - 


clockwise increase of 0^ 


. --/C ) decreases, and as 0 decrease&--counter ■ 

2 V P^w V 

clockwise increase of 0 -r 


3 ’/-fC \ increases. Then (C ) is a function 
2 ' p'w ' p'w 

of 0 for any particular piece of equipment. 


V V 

We can plot calculated (^p) 0 or (Vw vs 0^/tan p and 


experimental (K . ) vs 0 for each pump. We can then determine 
v' ^ 

(K . \ and (C \ at specific values of 0 or 0 /tan B for each pump. 
^ c, min'w ^ n'w ^ ^2 


p'w 
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The procedure for developing (K . ) vs ( C ) data are 

c, min w p w 

illustrated for three helical inducers on figure 5.4. Data similar to 
those shown in figure 5,4c are required for each value of ijj ^ of interest. 
Extensive compilations of such data may make it possible to derive 

V 

functional relationships between (K . ) and (C ) for certain 

c, min w p w 

classes of pumps. Interpolation or extrapolation of these data to new 
or different designs may then be possible. 

V 

If desired, (C ) and (K . ) can be expressed in terms 
' p w c, min w 

of the relative velocity, w , for the mean effective blade radius. 

m 

0 changes continuously with radial position along the blade leading 
edge for constant N, V , and 3 . Under these conditions 0 is minimum 

V ^ 

at the blade tip, ( is a maximum and thus cavitation should originate 

at the blade tip; therefore, the tip radius was selected as the pump charac- 
teristic dimension in this study. 

5.3 Estimation of Pump K . Using K . Data for Stationary- Bodies 

c, min c,min 

(Pump experiments are not required) 


This method of estimating K . is related to the K . vs 
^ ^ c, min c, min 

C approach outlined in subsection 5* 1 but is considerably more complex. 

p 

The following discussion was postponed until now so that the reader 
could benefit from the foregoing arguments and thus avoid confusion. 

On figure 4. 1 the NBS stationary-body data are adequately 
represented by 


(k . \ = 0. 30 C . (5-3) 

I c , min / CO p 

' V 

Due to the uncertainty in C for zero-caliber ogives Billet's [28] data 
were not used in determining the slope of this line. If the Billet 
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{ Kc, min) w 


y 

a) Calculated noncavitating (Cp)w 
( independent of 


Figure S. 


b) Experimental (Kc, min)w 

(f' = constant). 



c) Experimental (Kc,min)w and calculated (Cp)w 
( ^ = constant ). 


Development of (K 
( * - con stant). 


. ) 

min w 



data for a family of helical induce 




data are forced to fit this straight line, we obtain C = 1. 26 and 

P 


6 ) 


= 1. 20 for his zero-caliber ogives. Using the area correction 


factor, F , eq (5-3) can be replaced by the equivalent expression 
C 


2 _ 

K . 0.30 

C c, min 



(5-4) 


Equation (5-4) assumes that F is about the same for cavitating 
and noncavitating flows. This is a fair approximation for the relatively 
thin cavities developed on the stationary-bodies and for average cavity 
lengths associated with i^iin * P”iimp blade passages, the cavity 
thickness is frequently many times the blade thickness; therefore, the 
cavity shape can dominate F for cavitating pumps. Accordingly, 
must account for the presence of cavities when applied to pumps. 

In deriving the data represented by eq (5-4) it was noticed that 

2 — 

K . constant for all of the stationary-body data. Because 
c , mm 

K . from Billet is not in doubt, his data were included with the 
c, min 

NBS data to obtain 



c, min 


0. 374 . 


(5-5) 


An expression of this variety has great appeal because F inherently 
accounts for some of the detailed blade -plus -cavity geometry of each 
individual pump. This appeal exists even though the constant in eq (5-5) 
may be expected to vary for different classes of pumps and for different 
values of t|r^. Assuming the applicability of this expression to pumps. 



this computation, we will use the basic cascade configuration illustrated 
in figure 5. 5, This simple configuration is adequate for the inducers 
correlated herein but will require some modifying assumptions for 
application to the impellers. 
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Figure 5. 5 Cascade configuration illustrated for computation of F i 
cavitating inducers and impellers, ^ 






To evaluate F^, we must determine the blocked and unblocked 
areas of the blade passages as they appear to an observer traveling 
with a particle of fluid along the velocity vector w^ • The unblocked 
area is taken as 


( 


TT D. 


n 


Sin Y 



D 

t 



(5-6a) 


and the blocked area is described by 


where g = 



By definition 



(5-6b) 


(5-7) 


To evaluate 6 we use the expression for cavity shapes derived from our 

V 

ogive studies ( see discussion in section 4, 3, 3 of this report), 

0 75 

6 - 0.41 X • . (5-8a) 

V 


Converting 6 
obtain 


and X in eq (5~8a) from millimeters to centimeters, we 

0 75 

6 = 0. Z31 X , (5-8b) 

V 


where x cannot exceed 't/Z , 

Equation (5-8b) must be adjusted to account for the effect of the 
modified cascade solidity on cavity thickness. In our preliminary hydro- 
foil and ogive tests we observed that cavity thickness decreased almost 
linearly with increasing solidity (blockage) as it is normally defined in 
turbomachinery design. To adapt eq (5- 8b) to the pump data, we note 
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that the noncavitating blockage factor, F , for the stationary-bodies is 
analogous to t^/ S for the pumps. for the quarter- caliber ogives was 

0 . 095 , nearly the same value as for impeller B (see table 5.1). 

Thus, we will assume that eq (5-8b) is applicable to the prediction of 
cavity shapes within the blades of impeller B and will linearly correct 
this equation for solidity effects in the other equipment. The most 
meaningful noncavitating solidity parameter at our disposal is the modified 
cascade solidity at the tip, C . From table 4.4 we obtain C ^ = 1.683 
for impeller B and then eq (5- 8b) may be rewritten in the generalized form 


6 =0.231 

V 


/l.683\ , ,0.75 


(5-9) 


Equations (5-7) and (5-9) may now be used to evaluate F for 
all of the NASA-LeRG pumps. Any reasonable value of x may be 
inserted into eq (5-9)“-“the unblocked blade chord length, G^^, was used 
herein. All of the geometry and flow parameters necessary to evaluate 
these two expressions are contained in tables 4.4 and 5.1. Most of 
the geometry data given in tables 4,4 and 5. 1 are readily calculable 
according to the illustration given in figure 5.5; however, some of the 
impeller data required measurement of the equipment and/or reference 
to construction drawings. 

In table 5.1, the values of jg for the impellers are not the true mea- 
sured values but are calculated as shown on figure 5.5 and in the footnote 
of table 5.1. The calculated values were used after careful study of a 
NASA cavitation film [30]. This film, for a pump similar to impeller B, 
revealed that vapor always fills the space adjacent to the suction surface 
of the impeller blades (see figure 5.6). From this observation and for 
simplicity of analysis, we are effectively straightening the impeller 
vanes to conform to the illustration in figure 5.5. For these impellers 
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Table 










the cavity is assumed to have the additional thickness--measured from 

the theoretical straightened vane position- -as specified by eq (5-9). 

Again, this concept is illustrated in figure 5.6. 

Equations (5-5), (5-7), and (5-9) were applied to the NASA-LeRC 

impeller and inducer data and the results are presented in table 5, 1. 

The predicted values of K . are remarkably close to the experi- 

c, min •' ^ 

mental values obtained from our previous calculations. These predicted 
values of K . do not correspond to known values of il;^as explained 
in the following paragraph. For simplicity, inlet blade tip dimensions 
were used to evaluate the geometry parameters in this two-dimensional 
cascade flow illustration. In reality, the flow is much more complex 
(particularly for the impellers) and could be treated to account for varia- 
tion in blade height, curvature, thickness, etc. Also, mean effective 
dimensions might be formulated to produce better results in the simpli- 
fied analysis. Only a modest effort was expended in this direction be- 
cause so many possibilities exist; arithmetic mean values of blade thick- 
ness, modified cascade solidity, ig, etc., were tried and they produced 
much higher values of K . than those listed in table 5.1. Perhaps 
the initial choice of inlet blade tip dimensions was fortuitous because it 

is apparent that prediction of K . from eq (5-5) is highly sensitive 

c, min ® 

to choice of physical dimensions. More detailed analysis is needed to 

improve K . predictions by this method, 
c, min 

Possibly the major deficiency of this approach is that we have 

no idea what value of corresponds to our predicted min’ ’ 

we must determine what head loss is attributable to a cavity of half- 

length C . Obviously this head loss, and therefore will vary with 
ub _ IV 

each pump. Theoretically, — ► I as jr' — ►!. 0. Also, 

S p. S as i|f' ►!. 0 . To pursue this idea, we must establish 

c nc 
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relationships between K . and for each specific pump. At least 

c , min 

two potential methods exist for accomplishing this task- -both are 
described below. 

This blade -plus -cavity approach may be salvaged by performing 

idealized flow calculations to evaluate the cavitating head loss. These 

computations are performed at constant N and Q for a pump of known 

geometry. Cavity shape is predicted by eq (5-9) for an assumed cavity 

length and a corresponding value of niin calculated. Head loss is 

then computed by assuming that the cavity behaves like a solid body, i.e. , 

by assuming that the blades have the shape of blade-plus-cavity. This head 

loss is converted to for identification with its companion K . * 

^ c , min 

Such computations are not simple, but even rough estimates should 

be adequate to establish the needed relationship between predicted 

K . and Computerized solutions based on incremented cavity 

c , min ^ 

lengths may be terminated when the desired i|;^is reached. The cor- 
responding cavity length is then used to evaluate rnin previously 

described. Next, K . is used to evaluate pump NPSH requirements 

c, min 

using eqs (4-16) and the method described herein. 

More experimental pump data and analytical work are required 

to confirm or invalidate this approach. 

A second method of linking K . and ilr' requires experimentally 

® c , min 

verified^ vs tlr' data . Knowing the pump geometry and using eqs (5-9), 

V 

(5-7), and (5-5)>we can calculate 6 , F^, K , and for each 

V C c,min V _ 

assumed value of x; therefore, functional relationships between K 

c, min 

and may be derived for each pump. By relating^^ to we can 

establish the needed K . vs ilf' data. Experimental data, with intra- 

c , min 

blade flow visualization, is required to establish the ^ vs function 
for each class of pumps. Thus, this method of evaluating ^ has no 
appeal unless a simple and universal correlation between and ijf' 
is discovered. 
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In eq (4-21) we reasoned that'V^ « (l-i|r^) . We unsuccessfully 

attempted to relate K . to 'V' in our stationary-body tests. K 

c,min V r y . c,min 

was found to be independent for the venturi and hydrofoil and was 

weakly related for the ogives (K . ccV'"^*^^). From Billet’s 

c,mm v_Q ’ 

data [28] , it was deduced that K . <xA^ * , Substituting C (l-\lr^) 

c, min V o 2 \ T / 

for V , we may write 


or 



(5-lOa) 


K 


c, min 


(^c, min) ref 



(5«10b) 


where the body or p\imp geometry is held constant. To test this simple 
expression, values of z were computed for the impellers and inducers. 

These values were just as inconsistent as those for the stationary-bodies, 

^ — z 

Thus, the assumption that K . ccAf fails and eq (5- 10b) cannot 

c , min V 

be used to transfer K data from one to another, 

c, min 

The prospects of acquiring experimental verification of universal 

functions that relate V and are currently dim. Nevertheless, this 

method of coupling if' and K . is sufficiently- obvious to iustifv the 

c, min 

effort expended herein. 
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6. SUGGESTIONS FOR FUTURE WORK AND DATA COMPILATIONS 
This study has revealed a number of areas where future work 
is needed. Each of the tasks listed below will require considerable 
effort from the scientific community and pumping machinery specialists. 

1) Improve the accuracy and consistency of fluid thermophysical 
properties- -particularly for conventional fluids such as water, 
hydrocarbons, and the common refrigerants. 

2) Launch a significant effort to collect existing pump data (un- 
published and published), together with geometric details of 
pumps, for purposes of correlation and to test correlative- 
predictive procedures outlined in this report. 

3) Provide the type of data indicated in this report and provide 

the maximum in geometric detail for all future pump experimen- 
tation. Expend maximum effort to obtain data of the highest 
possible quality and precision. 

4) Evaluate the potential methods of estimating K . outlined 

^ ^ ® c,min 

herein and establish C - K . - 0 relationships with pump 

p c, min 

data compilations. 

V 

5) Improve and simplify, if possible, calculations. 

6) Perform additional stationary-body work because this study 
and previous work [3] show the relevance of such work and 

its direct applicability to the prediction of cavitating performance 
in rotating machinery. 

7) Finally, as a result of additional knowledge and experience, 
improve the existing correlative expressions. 
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7. CONCLUDING REMARKS 

The work performed in the course of this study is documented 
in the four volumes of this report series and may be succinctly sum- 
marized as follows: 

1) A simplified, more precise computation [23] was devised 
for the calculation of B-factors. This improved method 
is based upon isentropic vaporization of liquid and clearly 
demonstrates that precise thermophysical properties of 
fluids are required for improved future work. 

2) Convective heat transfer and two-phase mass flux limiting 
concepts were used [21] to extend and improve existing 
formulae for correlating developed cavitation data. These 
correlative expressions were extended to the prediction of 
cavitating performance of pumps in this final volume (IV), 

3) Extensive developed- cavity data were acquired for three 
different hydrodynamic body shapes using liquid hydrogen 
and liquid nitrogen test fluids. These data were obtained 
using two-dimensional and axisymmetric bodies and for both 
internal [20] and external flows. Size effects were obtained 
for the submerged axisymmetric bodies [22], 

4) The results of work performed in items (1) and (2) were used 
to correlate the data obtained in item (3). 

5) The liquid hydrogen pump data generated by NASA-LeRC 
were analyzed and correlated using the results obtained 
from items (1), (2), (3), and (4). Thus, the state-of-the-art 
of predicting cavitating performance of pumping equipment 
was significantly extended. 


73 



6) The relevance (to pumps) of analysis and experimentation 
with stationary bodies was strengthened; i.e., the 
correlative expressions derived for stationary bodies 
appear extendible to rotating equipment. 

7) The effort expended in items (1) through (6) revealed 
certain technical areas where additional analysis and 
experimentation are needed- -these areas are defined. 


74 



0 <! 


8. NOMENCLATURE 


adiabatic acoustic velocity of liquid, evaluated 
at the minimum cavity pressure 

adiabatic acoustic velocity of vapor, evaluated 
at the minimum cavity pressure 


A 

o 


B 


BELASH 


B 


t 


C 

n 


inlet area of test section (tunnel); may also be 
inlet area in line-mounted pump impellers or 
pump inducers 


ratio of vapor to liquid volume associated with 
the sustenance of a fixed cavity in a liquid 


[ .V /VJ 


B derived from isentropic flashing theory (Ref. [23]) 


BELASH 


chord length of inlet blade at hub, see table 4,4 

(n = 0, 1 . . . etc, ): constants or numerical 
coefficients in various algebraic expressions 


C 

P 


P 


(C ) 


p w 


noncavitating pressure coefficient 

minimum noncavitating pressure coefficient 

[ E ^ - h^/(v//2g^)l 

minimum noncavitating pressure coefficient 
based on inlet blade tip relative velocity 
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c 


t 



D 


D 


D 


D, 


D. 


h,ps 


D. 


h, ss 


D 


ml 


D 


m2 


D 


m3 


value of C in an infinite flow field fs F fC )] 

P C P 

specific heat of liquid at constant pressure, 
evaluated at the tunnel (pump) inlet (P and T ) 
modified cascade solidity at tip of inlet blade [ 
note that this solidity parameter is equivalent to the 
conventional definition of cascade solidity divided by 
cos 3 

chord length of inlet blade at tip, see table 4.4 

unblocked chord length of inlet blade at tip, see table 4.4 

characteristic dimension of stationary-body or 
pump equipment; taken as radius of contour for 
stationary-bodies and as inlet blade tip diameter 
for impellers and inducers 

discharge (overall) diameter of impeller 

hub diameter of inlet blade at the inlet 


hub diameter of primary splitter vanes on impellers 


hub diameter of secondary splitter vanes on 
impellers 


inlet blade mean effective diameter 
for impellers or inducers 


> 


impeller mean effective diameter 

l2 ;; Tz 


Z Z 

nD "l-nD “tZnD 

_t t,ps t , s s 


4n 


impel ler mean effective diameter 

_ 2 " 


2 2 

nD +nD +2nD 

ml m, ps m, ss 


4n 


I 2 2 

^t + ^ 
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D 

m, p s 


D 

m,ss 


D 

o 


D 


t 


D 

D 

F 


t,ps 
t, ss 
B 


F 


C 


Fr 


g 


c 


h 

c, min 
h 

o 


h 

V 


mean effe ctive diameter of primary splitter 

2 


vane s = 


[■Vi 


ps n,ps 


mean effective diameter of secondary splitter 


vanes 


[■V3 


SS 


+ D. 


h, SS 


tunnel inlet diameter for stationary-bodies or 
inlet duct diameter for impellers and inducers 

tip diameter of inlet blade at tbe inlet 

inlet tip diameter of primary splitter vanes on im- 
pellers 

inlet tip diameter of secondary splitter vanes on 
impellers 

area blockage factor, the fraction of cross sectional 
flow area occupied by the stationary-body or by the 
pump blade -plus -cavity 


area correction factor, the square of the ratio of 
blocked cross sectional flow area to unblocked (inlet) 
cross sectional flow areal 


Froude number 




r / 

n j 


conversion factor in Newton’s law of motion 
(gravitational acceleration) 


minimum static head in cavitated region corresponding 
to minimum absolute static pressure within the cavity 


tunnel or pump inlet static head corresponding to 
absolute static inlet pressure 

head corresponding to saturation (or vapor) pressure 

of the test liquid at the tunnel (pump) inlet temperature, 
T 

o 
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h 

X 


V 

h 


All 

V 


H 

AH 

c 


static head corresponding to absolute static pressure, 
measured on the stationary-body or on the suction 
surface of the pump inlet blade at distance x, down- 
stream of the minimum pressure point--for non- 
cavitating flow 

static head corresponding to the minimum absolute 
static pressure on the leading edge of the stationary- 
body or on the suction surface of the pump inlet blade, 

V 

computed from expression for C 

static head depression attributable to the vaporiza- 
tion of liquid that is required to sustain a developed 

cavity = h - h . 

^ V c,minj 

total head, corresponds to total pressure 

total head rise across cavitating pump impeller or 
inducer 


AH 

nc 


total head rise across noncavitating pump impeller 
or inducer 


k 


K 

c, min 


K 

c, min 


thermal conductivity of liquid, evaluated at tunnel 

inlet or at pump inlet (P and T ) 

o o 

developed cavitation parameter, based on minimum 
static cavity pressure 

r= /h - h \/(y “/2g \ 

I o c, min ] I o J 

arithmetic mean value of K . for a complete 

c , min ^ 

set of data points for a particular hydrodynamic 

body-fluid combination or for a particular pump 

at mean values of also used for arithmetic mean 

value of K .at specific operating conditions for 
c , min 

a particular pump 


78 




K 

V 

I 

MPE 

MTWO 


n 

n 

ps 

n 

ss 

N 

NPSH 

Pe 

PE 

P 

o 


value of K .in an infinite flow field 

[-c«'Z)] 

developed cavitation parameter, based on inlet 
blade tip relative velocity 

conventional developed cavitation parameter 

r , . , . 2 . -| 




length of cavities developed on stationary-bodies 
or on pump blades, used interchangeably with x 

arithmetic mean percent error in predicted values 
of NPSH- -see tables 4* 5 through 4. 7 


liquid phase velocity ratio | = V /V 


io = 






reference [21] 


number of inlet (main) blades on impeller or 
inducer 

number of primary splitter vanes on impeller 

number of secondary splitter vanes on impeller 

angular velocity of pump impeller or inducer 

2 

net positive suction head [= h +(V /2g )-h ] 

o o C V 


Peclet number = V t/a 


peripheral extent of blades, degrees of arc 


tunnel or pump absolute static inlet pressure 


79 


Pr 


Q 


Re 

Re 


D 

I 


>6 


s 

s 


fl 

f2 


S 

nc 

S 

c 



T 

o 


u 

t 


Prandtl number 




volumetric flow rate (capacity) of pump 


diameter Reynolds number - 


length Reynolds number 


z P V 
o o 

V t/ J 
o o 


peripheral length along the contour of a stationary- 
body as measured from the leading edge of the 
contour to the minimum pressure point on the 
contour 


specific entropy of saturated liquid at 

specific entropy of saturated liquid at minimum 
cavity pressure (temperature) 

specific entropy of saturated vapor at minimum 
cavity pressure (temperature) 

noncavitating suction specific speed of pump 

cavitating suction specific speed of pump 

the minimum distance (with allowance for blade 

thickness) between impeller or inducer blades as 

measured between lines constructed tangent to the 

blade pressure surfaces at the blade tips, see figures 

5 . 5 and 5 . 6 (tt D /n) sin 3 “ t, ] 

t b 

thickness of all blades on the impellers: tip 
thickness for inducer blades 

bulkstream static temperature of liquid entering the 
tunnel or pump 

pump inlet blade tip speed 
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V 


I 


V 

o 


characteristic liquid velocity component, normal 
to cavity liquid- vapor interface, see reference [21] 

bulkstream velocity of test liquid at inlet to tunnel 
or pump 


VHO 



V 


We^ 


w 


t 


X 


Greek 


experimental developed cavity data for the combined 
venturi, hydrofoil, and ogive stationary-bodies 

volume of liquid associated with the sustenance of 
a developed cavity 

volume of vapor associated with the sustenance of a 
developed cavity 

Z 

Weber number [ = p V ol 

o o 

liqui d velocity r elative to pump inlet blade at tip 



axial distance measured from minimum pressure 
point on stationary-bodies; also used to measure 
cavity length along impeller and inducer blade 
chords 


O' 


3 


Y 


thermal diffusivity of liquid, evaluated at tunnel 

(pump) inlet (P and T ) 
o o 

inlet (main) blade angle at tip--measured from 
circumferential direction, see table 4.4 and figure 5.3 

blade leading edge fairing angle [ = 3/2] 

-1 

tan 0 
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6 


V 

9 

2 


V 


thickness of the developed vaporous cavity- 
angle between blade suction surface and 
angle between blade pressure surface and w^ 

absolute viscosity of liquid^ evaluated at tunnel 

(pump) inlet (P and T ) 
o o 

kinematic viscosity of liquid, evaluated at tunnel 
(pump) inlet [ r. \i/ p^] 

density of liquid, evaluated at minimum cavity 
pressure 


P 

V 


a 


0 


i/ 

1 /' 


density of liquid, evaluated at tunnel (pump) inlet 

(P and T ) 
o o 

density of vapor, evaluated at minimum cavity 
pressure 

surface tension of liquid in contact with its vapor, 

evaluated at tunnel (pump) inlet (P and T ) 

o o 

flow coefficient Ff-V /u 1 
^ o t 



ratio of cavitating to noncavitating head rise 

coefficients [= AH /AH 1 
c nc 


Subscripts 

o = denotes tunnel (pump) inlet location 

ref = reference run (data point), or test conditions, to 

which a computation is being referenced when 
attempting to correlate pump cavitation performance 
via eqs(4-14)or eqs(4-l6) and when using 
eq (4-1 ) or eq (4-2) 
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Superscripts 

El 


E2 


E3 

E4 

E5 

E6 

P 


exponent on thermal diffusivity ratio in eq (4-1) 
and eq (4-2) 

exponent on tunnel inlet velocity ratio in eq (4-1) 
and also used as an exponent on the MTWO ratio 
in eq (4-2) 

exponent on cavity length ratio in eq (4-1) and 
eq (4-2) 

exponent on kinematic viscosity ratio in eq (4-1) 
and eq (4-2) 

exponent on surface tension ratio in eq (4-1) and 
eq (4-2) 

exponent on characteristic dimension ratio in 
eq (4-1) and eq (4-2) 

exponent in algebraic expression for cavity shape 

(6 = C xP) 

V o 

exponent in algebraic expression, defined by 
eqs (5-10) 
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APPENDIX A: Examples of NPSH Prediction Calculations t 


Examples of the calculational procedure for predicting NPSH 
requirements of cavitating impellers and inducers are given. For 
convenience, eqs (4-14) and (4-16) are repeated below: 



(4-16) 


Lacking cavity length data, the ratios in eqs (4-14) and (4-16) are 

taken as unity-- see discussion in section 4.6. 

The definitions of BFLASH and MTWO are also given as 


B ^ BFLASH =1 


I 


V 


s - s 
fl f2 


(A-1) 


v2 


fl 


V 

MTWO = — 
a 



0. 5 


t Numerical examples in this Appendix were supplied from the 
computer computations performed by William R. Parrish. 


(A-2) 
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Computational Procedure 

Predicting pump cavitation performance using eqs (4-14) and 

(4-16) requires two steps. First, from existing pump cavitation data, 

calculate reference parameters. Then, it is simple to predict a NPSH 

at different conditions as long as is held constant. 

Reference Parameter Calculation ; For example, say that we 

want to use eq (4-16) to predict NPSH. It is assumed that K . is 
^ ^ c, min 

some function of 0 and \|;' for each impeller and inducer. Two experi- 
mental data points are required to calculate reference parameters--at 
least one of these data points must exhibit measurable 'head depression'. 
Ah . This criterion virtually excludes the use of cold water as test 

V 

fluid for more than one test point. For each pump we may choose two 
data points at the same values of 0 and \lr^--under these conditions 
eqs (4-16) simplify to 



NPSH -h Ah 
(NPSH + 


(A-3a) 


and. 



^ MTWO 

MTWO , 1 
ref j 

< / 


(A-3b) 


There are six unknowns: Ah , (Ah ) . , B, B , MTWO, and MTWO , 

V V ref ref ref 

The relationships between Ah and B, and between B and MTWO are 

V 

nonlinear and vary with fluid, fluid inlet temperature and inlet velocity. 
However, these nonlinear relationships can be linearly interpolated 
over sufficiently short increments of B, 
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The computational procedure is as follows: 

1) Guess a value of determine a corresponding 

value of from an appropriate table or graph [23], 

Equation (A-1) may be used to compute directly using 

a temperature- entropy diagram for the appropriate fluid. 

2) Calculate Ah from eq (A- 3a). 

V 

3) By iteration, compute the correct values of B and MTWO 


satisfying eq (A- 3b). 

4) Interpolate [23] the value of Ah corresponding to the B 
found in step 3. 


5) Compare the values of Ah calculated in steps 2 and 4. 

; V 

If they are equal the correct value of (Ah ) , has been 

V ref 

» chosen, otherwise choose another value of (Ah ) and 

I V ref 

I return to step 1 . 


The difference in the calculated values of Ah in steps 2 and 4 

V 

is nearly linear in (Ah ) Therefore, linear extrapolation or 

V ref 

interpolation to zero difference is a good guide to choosing the next 


value of (Ah ) Similarly, in step 3 the difference between the left 

V ref 

and right sides of eq (A-3b) vs B is nearly linear; therefore, improved 

guesses for B can be made in the same way as described for (Ah ) 

v ref 

Predicting Pump Performance : With reference parameters 

known, predicting pump cavitation performance for another fluid, 
fluid temperature, rotative speed, flow condition, or pump geometry 
at the same value of il; ^ is straight-forward. B is computed from the 
second equation of (4-16) by iteration as described in step 3. The 
corresponding value of Zh is then determined as in step 4, and sub- 
stituted into the first equation of (4-1 6). Assuming that K . and 

c, min 

(K . ) ^ are known, the first equation of (4-16) is then used to 

c, min ref 
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solve for the desired (predicted) NPSH. An example illustrating the 

above procedure is included at the end of this appendix. 

Data Reduction; Determining which equation set works best 

involves the following procedure: 

1) The pump data were divided into sets by holding pump 

geometry, speed (inducer data only), 0, and constant. 

The head depressions for two data points are then calculated 

using the procedure previously outlined in this appendix. 

Using these calculated values of head depression, the 

corresponding values of K . and S are calculated from their 

c,min c 

definitions. Repeat this procedure for all possible pairs 

of data points within the data set. Arithmetically average 

the calculated values of K . and S and assign the 

c , min c 


average values to each data point in the data set. K 


c, min 


vs 0 and/or S vs 0 have thus been determined. 

Sequentially use all possible combinations of pairs of data 
points to determine reference data for each pump geometry 
and value of \|f^ Use each pair of reference data points to 
predict NPSH for all other data points for all other pump 
geometries at the chosen value of Finally, tabulate the 

average error in predicting NPSH for each geometry and 


This procedure is applied to each set of predictive equations to 
find the set which gives the lowest mean percent error (MPE) in 
predicted NPSH for each pump geometry, e. g. , the equation set giving 
the lowest MPE when predicting NPSH values for the impellers and 
inducer B using inducer A for reference data. 
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Numerical Example 

As an example of how to use the MTWO equations, we will predict 

the required NPSH for inducer B using reference data from impeller A. 

The equations are based on a constant and the inducer data 

are at i|;' 0. 7--so choose data points 7 and 8 (table 4. 3) to obtain 

reference parameters for predictive purposes. 

Calculation of Reference Data: To use the MTWO equations, 

reference values for B, MTWO, and Ah are required. If K is 

V c, min 

known (Ah ) can be calculated from the definition of K : however. 

V ref c, min 

it is assumed that is unknown for this part of the example. By 

choosing data points 7 and 8 (0 = constant) we may use eqs (A- 3) to 
evaluate the reference parameters. These simplifications do not alter 
the method of solution because the remaining terms in eqs (4-16) produce 
ratios of unity under the prescribed conditions. By calling data point 7 
the reference point, eq (A- 3a) becomes 


+ /97Ann\2 

(tw) <a-4, 

As a first guess, assume ~ 0 m of liquid hydrogen. 

From (A-4) Ah is 32, 92 m. This value will be compared with Ah 

V V 

found from the B-factor equation- -the B for the test point must be 

found from eq (A- 3b). Calculate B for the reference data point 

using eq (A-1), For (Ah ) = 30, 0 m the BFLASH is 0, 8185, With 

V ref 

B known, MTWO as calculated from eq (A-2) is 0,3817 . Then eq (A-3b) 
can be written as 

g = / MT^\ 

0. 8185 \ 0. 3817 j 
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In principle, eq (A-5) can be uniquely solved using the definition 
of MTWO; however, it is easier to guess a value of B, calculate MTWO 
and test this value of MTWO in eq (A-5). As a first guess, assume 

B = 0. 83--the corresponding value of MTWO is 0. 426, Putting this 

( 0 426o\®* 

o' 38 ~ i 7 J ^ • Comparing the 

assumed and calculated values, the initial guess of B was too small. 

/o 457 r\ 

Now guess B = 0. 942, then MTWO = 0. 4578 and B = 0. 8185f ^ 3817 j 
= 0, 898, This time the assumed value is too large; try B = 0, 878 . 

/O 4392\0. 51 

This corresponds to MTWO 0. 4392 and B = 0. 8185 [ * j 


= 0. 879; this calculated value of B is close enough to the assumed value 
of 0. 878. 

This value of B implies a head depression, via eq (A-l),of 31, 76 m 

which is smaller than the 32. 92 m calculated from eq (A-4). Therefore, 

the assumed value of (Ah ) ^ was too high; now try (Ah ) , 20. 0 m. 

V ref V ref 

From eq (A-4) Ah is 20. 815 . B-factor for the reference point is 

^ V 

0. 5035 and (MTWO) is 0, 2893. Repeating the same iterative pro- 

ref 

cedure described above, B-factor for the test point is 0. 5405, MTWO is 
0. 3325 and Ah is 21. 271 m. This time the head depression calculated 

V 

from eq (A-1) is larger than the value calculated from eq (A-4). 

The next guess of Ah . is found, by linearly interpolating to 

V ref 

zero difference in calculated head depression, to be 22. 68 m. From 
eq (A-4) Ah is 24. 00 m. The B-factor for the reference point is 

V 

0, 5830 and (MTWO) is 0, 3140, Solving eq (A-3b) by iteration gives 

ref 

a B-factor for the test point of 0. 6260 which corresponds to an MTWO 

of 0. 3611 and Ah of 24. 11 m. Therefore, the correct value of (Ah ) 

V V ref 

is 22. 68 m, (B) , = 0. 5830 and (MTWO) , = 0. 3140. 

ref ref 

To obtain reference values for eqs (4-14), the procedure is the 
same as above except that no iterative procedure is required to cal- 
culate the test point B-factor. 
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Prediction of NPSH at Different Rotative Speed, Temperature, 

Flow Coefficient, and Pump Geometry ; Using the reference data, 

predict the required NPSH for inducer data point 54. It is assumed 

that K . is known (or can be estimated) for the reference point and 
c, min 

test point. Referring to table 4. 3,eqs (4-16) become, for this problem, 


1. 0 + 3. 038\ /O. no (30000)12. 65 


1. 0 + 1. 402 /Vo. 225 (25000) 6. 78 


NPSH + Ah 


(21. 3 + 22. 68) 


or 


NPSH + Ah = 88. 58 m; 


(A- 6a) 


B / MTW0 \°* / 12. 65 \°* 

0. 5830 " \^0. 3140 f 6. 78/ 

or 

B = 1. 6389 (MTWO)°* . (A-6b) 

First eq (A-6b) must be solved by iteration as previously described. 

Guess B to be 1. 31; this requires that MTWO = 0, 7229 for this test point. 

0 51 

Then B = 1. 6389 (• 7229) * = 1. 389, which means the guessed value 

is too low. Try B = 1, 458, which gives an MTWO = 0. 7708 and 

0 51 

B = 1. 6389 (0. 7708) * - 1. 435; this value implies that the guessed 

value is slightly too high. Now try B = 1. 424 and correspondingly 
MTWO = 0. 7597, then B = 1. 6389 (0. 7597)°* = 1. 424; this is the correct 

value and corresponds to a Ah of 21. 15 m. From eq (A-6a) NPSH 

V 

= 88. 58 - 21. 15 = 67. 43 m. From table 4. 3 the experimental NPSH 
is 67. 0 m. 

Again, if eqs (4-14) are used, the procedure is identical to the 
one just described except that no iterative scheme is required to 
calculate the B-factor. 
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The procedure described and illustrated herein may be applied 

to any pump geometry, fluid, fluid temperature, rotative speed, flow 

conditions, etc. , as long as ijr^ is held constant and appropriate values 

of K . (or S ) can be estimated. If pump geometry, and 0 

c, min c 

are held constant, K . (or S ) may be assumed constant and 

c , min c 

eliminated from the correlative- predictive expressions, i. e. , in 
eqs (4-14) and (4-16). 
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